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This  research  effort  addresses  the  fundamentals  of  noise  in  electron  beams 
in  crossed  field  microwave  tubes.  Recent  evolvement  of  new  applications  for 
DOD  requires  tubes  with  better  noise  characteristics  in  addition  to  the 
higher  bandwidth  and  gain  characteristics  of  crossed  field  devices.  Much 
work  was  expended  in  this  type  of  research  up  to  1969  but  thereafter  concepts 
have  remained  neglected  and  unexplored  until  1977  when  efforts  were  renewed. 
The  approach  her*  is  to  re-examine  the  analyses  using  simple  models,  to - _ — 
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resolve  various  unknowns  and  gaps  In  knowledge  and  to  investigate  original 
unexplored  avenues  with  which  we  can  gain  more  Insight  into  the  mechanisms 
of  fluctuation  and  modulation  in  electron  beams  travelling  in  crossed 
electric  and  magnetic  fields. 

In  the  first  year's  effort,  expressions  were  derived  based  on  kinetic 
theory  for  the  seven  moments  over  the  energy  distribution  function  in  a 
plasms  thermionic  diode  with  crossed  electric  and  magnetic  fields,  both  for 
the  space-charge  limited  and  magnetic-field  limited  situations.  The 
expressions  are  given  as  integrals  over  the  distribution  function  at  the 
cathode  surface,  whose  form  is  left  arbitrary.  Analytic  results,  valid 
in  the  region  of  the  potential  minimum,  are  provided  in  the  space-charge 
limit  when  the  potential  profile  is  assumed  to  be  parabolic  in  the  vicinity 
of  its  minimum.  The  seven  moments  are  the  density  profile,  the  two  current 
flux  components,  the  two  mean  velocity  components  and  the  two  mean  square 
velocity  components.  Prom  the  latter  four  moments,  the  velocity  noise  is 
determined.  At  the  potential  minimum,  the  velocity  noise  is  minimized  when 
ujp,,2/**),.2  ■  2  where (Up^Vs  the  plasma  frequency  at  the  potential  minimum 
and  b>c  is  the  cyclotron,  frequency,  and  operation  under  this  condition  is 
recommended.  It  is  also  shown  how  Lindsay's,  Ho  and  Van  Duzer's  and  the 
present  analyses  are  relAted  to  each  other. 

In  the  current  year's  effort,  we  concentrated  on  applying  the  moment 
macroscopic  equations  to  Investigate  the  propagation  of  fluctuations  in 
crossed  fields.  For  the  electron  gun  region,  the  differential  equations 
were  derived  for  accelerate^  beams  being  bent  by  the  magnetic  field.  As 
a  start  in  computational  analysis,  we  concentrated  our  attention  on  the 
interaction  regions.  Here  the  nature  of  straight  electron  beams  with 
finite  width  and  transverse  velocity  gradient  was  explored  using  the  model 
of  Buneman  et  al.  Curves  wer^  generated  of  dispersion  relations  using  a 
numerical  eigenvalue  scheme  and  compared  with  Buneman 's  analytical  results. 
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APPENDIX  I:  Reprint  of  Article 


The  contractor  shall  furnish  scientific  effort  needed  to 


conduct  research  directed  towards  the  following  areas 


(a)  An  electronic  equation  Involving  differential  equations  will  be 


derived  for  the  aodulatlon  quantities  In  a  linear  crossed-f ield 


M-type  device,  which  will  include  the  dc  variations  In  the  beam 


(b)  The  above  relation  will  be  applied  to  investigate  the  growth  of 


noise  in  the  transition  from  the  gun  region  to  the  slow-wave 


structure  region,  including  dc  space-charge  effects  near  the 


cathode 


(c)  The  relation  will  also  be  applied  in  the  drift  region  between 


parallel  plates  and  the  related  diocotron  instability  will 


be  investigated 


(d)  Parameters  in  the  theory  will  be  varied  with  the  aim  of  reducing 


the  noise  content.  Possibilities  include  smooth  and  abrupt 


velocity  variations  and  various  cathode  lengths  and 


(e)  Kinetic  theory,  using  the  electron  velocity  distribution  function 


and  the  Vlasov  equation  ,  which  includes  the  thermal  spread  of  the 


will  be  investigated  for  the  dc  beam  parameters  and  then  for 


the  above  four  tasks  on  the  modulation  quantities  in  the  beam 


(f)  A  hybrid  approach  to  investigate  noise  will  be  carried  out  in 


which  the  kinetic  theory  approach  will  be  applied  in  the  region 


between  the  cathode  and  potential  minimum  leading  to  noise  boundary 


condition*  at  the  potential  minimum  while  the  electronic 
equation  will  be  applied  In  the  succeeding  region. 

(g)  If  the  Indications  justify  the  following,  theoretical 

Investigations  will  be  performed  on  new  cathode  and  gun 
structures,  attempting  to  eliminate  the  above  mentioned 
problems.  These  techniques  Include  velocity  jump 
transformers,  shortening  the  cathode  length,  tilting  the 
cathode  with  respect  to  the  magnetic  field,  cold  cathode 
operation  and  artificial  means  of  controlling  or 
eliminating  the  potential  minimum  in  front  of  the  cathode. 

1.2  Work  Completed 

During  the  first  year  we  have  investigated  parts  of  tasks  d,  e,  and  f , 
mainly  part  f.  In  the  second  year,  our  effort  was  directed  towards  casks  a, 
b,  c,  particularly  c.  In  the  future  work,  we  propose  to  concentrate  on  tasks 
b,  d,  f,  and  g.  The  aacrospic  moment  approach  will  be  pursued  Intensively, 
since  it  is  simpler  to  apply  for  analytical  results.  Generally,  the 
objective  is  to  investigate  the  growth  of  noise  in  crossed  fields  as  the 
beam  propagates  through  various  structures,  such  as  the  gun  region,  drift 


region  or  circuit  region. 


2-0 


INTRODUCTION 


In  recant  years  microwave  tubes,  especially  klystrons  and  magnetrons, 

have  been  made  more  efficient  and  lighter^glving  solid  state  sources  strong 

competition.  Magnetrons  have  been  very  flexible  in  adapting  tube  designs  to 

2 

systems  requirements  .  For  example,  Varlan  Associates  now  has  a  1  MW  S-band 
tube  which  can  be  replaced  In  a  fixed  cavity  structure.  High  gain,  to  30  dB, 
has  been  achieved  in  CFA  (crossed- field  amplifier)  devices  by  Northrop  Corporation, 
while  bandwidth  is  over  half  an  octave  at  3  Kw.  Increased  gain  means  that  overall 
system  costs  and  complexity  can  be  reduced  by  using  fewer  driver  stages.  Tube 

lifetime  has  also  been  extended,  being  as  high  as  8-10,000  hours  for  distributed 

3  4 

emission  tubes  .  In  future  years,  reliability  of  tube  systems  will  be  increased 

and  the  frequency  range  will  be  extended  into  the  mm  range  as  gyrotron  tubes 

are  developed^.  .This  latter  tube  is  projected  to  yield  power  output  in  the 

multikilowatt  range\  and  may  be  applicable  to  solar-power  satellite-to-ground 

power  transmission^  as  well  as  energy  injection  into  tokamaks  for  fusion^. 

While  engineers  extend  lifetime,  power  and  frequency  range  of  microwave 

cools,  fundamental  research  into  the  nature  of  the  cathode  surfaces,  potential 

distributions  and  electron  beams  has  been  continuing  apace.  While  magnetron-type 

electron  guns  are  used  in  CFA  devices,  they  are  also  applicable  to  the  generation 

of  the  hollow,  cylindrical  beams  used  in  gyrotrons*.  These  beams  are  expected  to 

allow  40K  conversion  efficiency  because  of  optimized  interaction  with  the  RJF 

field.  Electron  gun  design  and  beam  voltage  are  therefore  critical  parameters. 

Beam  trajectories  have  been  modelled  on  computers,  both  for  gyrotrons^  (IEEE  Trans-ED. 

g 

issue  June  1977)  and  for  CFA  devices  .  With  such  codes  it  is  possible  to  alter 
beam-shaping  electrodes  and  study  the  effect  on  electron  beam  shape. 

Trajectory  analysis  necessarily  involves  a  close  knowledge  of  the  potential 

9 

distribution  within  the  gun  region  of  the  given  device.  Kino  (1960)  designed 
a  "short  gun"  with  a  cathode  short  relative  to  the  cycloid  length,  and  a  "long  gun" 


for  increased  power.  Masoarl  and  Rowe  (1965)  used  a  technique  developed 


and  Kirstein  to  calculate  the  electrode  shapes  and  potentials 


required  to  produce  a  laminar,  or  Brlllouin,  electron  beam  in  the  short  Kino 


gun.  More  recently.  True  (1978)  presented  numerical  analyses  of  potentials 


within  Kino  guns,  using  a  deformable  mesh.  Such  analyses,  although  in  two 


dimensions,  are  realistic  approximations  to  the  actual  conditions  within  an 


Fontana,  MacGregor  and  Rowe  (1979)  have  developed  an  improved  two 


dimensional  computer  simulation  of  a  Kino  gun,  using  5-10,000  representative 


electrons.  The  potential  distribution  Is  that  of  True 


and  Includes 


secondary  emission  from  the  cathode.  Two  methods  are  used  to  advance  the 


trajectories:  in  the  static  calculation,  convergent  solutions  are  found  between 


successive  solutions  of  Poisson's  equation,  while  in  the  time-dependent  calculation 


no  such  steady  state  Is  sought.  They  analyze  several  types  of  electron  guns 


Including  the  short  and  long  Kino  guns.  Convergent  solutions  were  found  readily 
for  the  short  Kino  gun  but  for  the  long  gun  returning  cyclolding  electrons  rendered 


convergence  unattainable.  For  this  latter  case,  the  time-dependent  model  proved 


successful 


Since  1964  much  empirical  work  has  been  done  on  reducing  the  noise  level  in 


the  long  Kino  gun.  In  view  of  their  light  weight,  high  power  and  wide  bandwidth 


reducing  the  noise  level  in  Injected  beam  devices  is  expected  to  go  a  long  way 
toward  improving  their  scope  of  use.  To  this  end,  Fontana  et  al  have  been 


developing  their  two-dimensional  model,  while  Kooyers  and  Shaw  (1979)  have 


produced  a  three-dimensional  model  using  a  slightly  simplified  geometry.  They 


have  developed  a  Lagranglan  particle  simulation  code  capable  of  modelling  up  to 


20,000  electron  charges.  The  electrons  are  emitted  from  the  cathode  with 


random  velocity,  number  and  position,  and  from  their  subsequent  positions 


and  velocities,  currant  and  valodty  fluctuations  (noise)  are  calculated.  They 
define  a  "space  charge  smoothing  factor"  T,  as  ratio  of  root  mean  square 
fluctuation  at  gun  exit  to  rms  fluctuation  at  the  cathode;  T  is  found  to  be 
'v  1.75  in  the  temperature  limited  case,  and  ^  0.2  In  the  space-charge  limited 

case. 

The  ultimate  source  of  the  noise  in  a  CFA  is  the  randomness  in  the 
emission  current  density  leaving  the  cathode  surface^.  Trajectory  modelling  is 
one  way  of  mapping  the  development  of  noise  from  cathode  to  gun  exit.  However, 
the  mechanisms  of  fluctuations  in  electron  beams  can  be  solved  beat  by  a 
consideration  of  the  motion  of  particles  in  cross-fields  near  the  potential 
minlmtmi  using  kinetic  theory.  On  the  other  hand,  the  propagation  of  noise  can 
be  most  readily  investigated  by  a  consideration  of  the  field  equations  associated 
with  the  beam  together  with  the  equation  of  motion. 

In  the  first  year's  work, we  successfully  analyzed  using  kinetic  theory 
Che  velocity  noise  in  crossed  fields  near  the  potential  minimum.  We  derived 
a  basic  condition  on  density  and  magnetic  field  to  minimize  the  mean  square 
deviation  in  velocity.  This  work  is  reported  in  Section  3. 

In  the  recent  year's  work, we  undertook  to  investigate  the  full 
equation  approach,  reported  here  in  Section  4. 


VELOCITY  NOISE  IN  CROSSED  FIELDS 


A  detailed  description  of  this  phase  of  the  vork  is  provided  in  the 
published  reprint  reproduced  in  Appendix  I.  Here  ve  briefly  summarize  the 


results 


He  investigate  both  the  space-charge  United  case  and  the  magnetic 


field  limited  case  in  a  thermionic  diode  with  crossed  electric  field  along  the 


x-direction  and  magnetic  field  along  the  z-dlrectlon.  By  a  careful  systematic 


analysis  based  on  an  assumed  parabolic  spatial  variation  in  the  potential 


about  the  space-charge  minimum,  ve  clarify  several  aspects  on  the  operation 


of  a  space-charge  limited  diode  with  a  transverse  magnetic  field  and  ve  obtain 


simple  but  Important  results  for  the  noise  in  the  velocity  components 


In  the  first  part  of  the  analysis,  ve  illuminate  Lindsay's 


(1960 


1964)  considerations  on  the  velocity  regions  over  vhlch  ve  have  to  integrate  the 


velocity  distribution  function  to  derive  the  moments  or  integrals  over  this 


distribution.  Let  subscript  zero  refer  to  the  initial  velocities  v  and  v  at 


the  cathode  upon  emission.  Without  subscript  zero,  the  velocities  (v  and  v  ) 


refer  to  those  at  position  X  vlthin  the  diode.  *  The  equivalence  betveen  the 


v  )  and  (v  ,  v  )  planes  as  far  as  the  regions  and  limits  of  integration 


are  identified  and  mapped  both  in  the  space-charge  and  magnetic  field  limited 


situations 


For  the  space-charge  limited  case,  Lindsay  (1964)  hovever  did  not  make 


much  headvay  since  his  integrals  could  not  be  manipulated  before  the  space 
potential  distribution  vas  known.  Working  in  the  (v  ,  w  )  plane,  ve  Illustrate 


a  number  of  new  features,  upon  assuming  the  potential  to  have  a  parabolic  spatial 


variation  about  the  potential  minimum.  (We  also  justify  this  assumption  by 


generalizing  somewhat  the  results) 


In  relation  to  the  emission  current  flux,  J  ,  perpendicular  to  the 


electrodes,  we  prove  that 


a)  our  formual  reduces  to  Lindsay's  (1964)  if  one  inserts  our 

assumed  parabolic  potential  variation  into  Lindsay's  equation, 

19 

b)  the  result  in  Ho  and  Van  Duzer  (1968) ,  also  for  a  parabolic  spatial 


variation  in  potential,  agrees  with  our  result  if  we  allow  the 


anode  to  receed  to  infinity,  as  they  do, 

c)  and  thereby  we  have  clarified  the  relationship  between  Lindsay's 
and  Van  Duzer' s  analyses. 


General  relations  are  derived  for  seven  moments  in  terms  of  integrals 


over  the  two-dimensional  electron  velocity  distribution  function,  both  for  the 


space-charge  limited  and  magnetic  field  limited  situations.  The  distribution 


function  refers  to  its  value  at  the  cathode  but  its  form  is  left  arbitrary 


except  for  the  assumption  that  it  is  a  function  of  the  magnitude  of  w 


density,  the  mean  velocities,  v  and  v  ,  and  the  mean  square  velocities,  v  and 


For  the  space-charge  limited  situation,  we  effect  major  simplifications 


upon  adopting  all  of  Ho  and  Van  Duzer's  (1968)  assumptions.  That  is,  besides 


using  a  parabolic  potential,  the  reasonable  assumptions  are  also  made  that  the 


potential  minimum  is  very  close  to  the  cathode  and  very  far  away  from  the  anode 


With  these  assumptions,  the  relations  for  the  moments  can  be  expressed  in  terms 


of  simple  integrals  over  the  three-dimensional  distribution  function,  which  can  be 


readily  integrated  for  a  Maxwellian  form 


The  cutoff  characteristics  of  J  versus  magnetic  field,  B,  differs 


markedly  between  the  magnetic  field  limited  case  and  the  space-charge  limited  case 


In  the  former  case,  J  depends  on  magnetic  field  and  anode  conditions  (position 


potential)  and  it  decreases  abruptly  above  a  critical  magnetic  field  given  by  the 


the  important  space-charge  parameter  uy  v.  —  -n -  m  m 

For  a  parabolic  potential,  6  Is  given  by  where  ^  1.  the  cyclotron 

frequency ,  u  Is  the  plasma  frequency  and  subscript  m  means  that  the  quantity 
is  to  be  evaluated  at  the  potential  minimum.  Space-charge  limited  conditions  are 
characterised  by  6>1.  whereas  magnetic  field  limited  conditions  occur  when  8  <1. 
when  8>1.  the  current  flu*  Jx  decreases  slowly  as  J.KS  -  D/B^e™  as  the 
magnetic  field  (B-S'*)  Is  Increased.  Also  Its  magnitude  Is  down  by  a  factor 
."m  from  the  temperature  limited  J.  value,  similar  to  space-charge  conditions 


gives  J  -  JgP  •  The  108  8um*  wx  y  s  » 

in  the  absence  of  a  magnetic  field.  This  dependence  for  Jy  Is  again  quite  different 
from  calculated  results  by  Lindsay21(1962)  and  Lindsay  and  Goodwll22(1965)  In  the 
magnetic  field  limited  situation. 

From  the  relationships  for  the  velocity  moments,  and  vx>y‘.  »«  obtain 

the  mean  squat,  daviations,  *"<*  ln  !'*rtlcul*t  th“lt  V,1““  “ 

potential  minimum.  The  noise  velocity  components  are  the  mean  square  fluctuations 
of  the  velocity  components  which  are  created  by  shot  noise  ln  the  emission.  Using 
Rack's  (19 38) 2 3  theory,  we  relate  the  velocity  noise  components,  \vx>y  X°  the 
_ _ w.ineicv  deviations .  The  results  are  illustrated  here  ln  Figure  1.1. 


QUANTITY  SHOWN  NORMALIZED  TO  leAfkTexp(-rVn)/(mJ*A#®)l 


Figure  3.1 :  Variations  of  Moan  Square  Deviations,  x  and  y  Components  and  Sum, 
Versus  S  ®t  the  Potential  Minimum. 


In  th«  absence  of  a  magnetic  field 


the  familiar  multiplicative  factor 


<4  —  ir>  -  0.8584,  (See  Smullin  and  Haua 


preaence  of  a  magnetic  field.  Van  Outer  (1963)  uses  the  magnetic  field 
independent  factors  (4  -  it)  and  2  respectively  for  (v  2>  and  <v  2>.  Our  results 


are  much  different.  Both  components  depend  on  magnetic  field  through  the  8 
factor.  The(v  component  is  large  and  uegatlve  for  low  8  values  near  one.  Its 


value,  normalised  as  shown  in  the  graph,  approach  4  -it  as  6  becomes  large.  On  the 
other  hand,  the  normalized^2)  starts  off  when  8  ■  1  at  4  -  it  and  decreases  to 
large  negative  values  for  large  8.  Considering  the  total  velocity  noise  content, 
namely  (v^2)  +(v^2),  the  graph  shows  that  its  magnitude  has  a  minimum  at  8  ■  2  where 
its  numerical  factor  is  2t/i(5-2it)  •  -3.629.  Although  this  is  4.23  times  worse 


than  the  zero-magnetic  field  value,  it  represents  the  least  velocity  noise  that  is 


obtainable  at  the  potential  minimum  position.  Thus  for  optimum  operation,  one 


should  allow  the  magnetic  field  to  satifsy  8  ■  w 


velocity  noise  is  also  proportional  to  the  cathode  temperature.  If  the  tube  can 


be  operated  wl.h  no  cathode  heater,  this  broadband  noise  will  be  greatly  reduced 
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This  probably  explains  the  experiments  of  Bespalov  et  al  *  (1967a, b). 


The  main  limitation  in  the  analysis  is  the  assumption  that  the  magnetic 


diode  is  either  space-charge  or  magnetic  field  limited.  In  most  cases  at  high 


current  emissions,  the  diode  is  in  fact  space-charge  limited  near  the  cathode  and 


potential  minimum  regions,  but  the  situation  changes  over  to  the  magnetic  field 


limited  case  in  most  of  the  diode  far  beyond  the  potential  minimum.  Our  result 


on  the  electron  density  profile  indicates  that  such  a  change-over  does  occur 


Our  initial  first  approximation  of  a  constant  density  profile,  inherent  in  a 


parabolic  spatial  potential,  is  found  to  be  Incorrect  in  the  second  approximation 
The  electron  density  is  found  to  decrease  rapidly  from  the  cathode  position  and. 


-Il¬ 


ia  the  second  approximation,  its  profile  is  similar  to  the  non-magnetlc  field 

cese.  If  we  define  sn  "effective"  6  -  uj  2  /ui  2 ,  we  note  thst  it  is  spetielly 

P  c 

dependent,  decreesing  from  6  >  1  to  8  <  1  at  positions  much  beyond  the 
potentiel  minimum.  An  anslyticel  solution,  incorporating  both  space-charge 
end  magnetic  field  United  situations,  is  non-existent,  but  an  iterative  computer 
solution  is  feasible  to  perform,  although  it  has  not  yet  been  done  by  us. 
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4.0  THE  FIELD  EOPATION  APPROACH 

4.1  Introduction 

Several  problems  can  be  solved  only  by  a  consideration  of  the  field 

equations  associated  with  the  beam,  together  with  the  equations  of  motion. 

This  field  equation  approach  has  been  undertaken  during  the  second  year’s  effort. 

The  noise  is  propagated  and  amplified  under  the  Influence  of  a  plasma  instability 

known  as  the  diocotron,  or  slipplng-stream  instability.  In  an  electron  beam 

drifting  in  crossed  electric  and  magnetic  fields,  a  growing  wave  exists  as  a 

results  of  the  velocity  slip  between  the  bottom  and  top  edges  of  the  beam. 

This  has  been  investigated  in  Idealized  situations  by  several  investigators, 

28 

notably  MacFarlane  and  Hav  (1950)  ,  who  treated  the  Brillouin  beam  (with  each 

electron  following  a  straight  path  and  b)  »  u  ,  where  u  is  the  gradient  of 

SC  s 

29 

velocity  across  the  beam  and  is  the  cyclotron  frequency) .  Gould  (1957) 

30 

treated  thin  beams  with  fa)  <<w  ,  Knauer  (1966)  analyzed  beams  of  finite  width 

sc 

31 

and  fa>  «  w  ,  and  Buneman,  Levy  and  Linson  (1966)  treated  finite  beams  with 
s  c 

arbitrary  fa)  :  o>  ,  ratio.  To  date,  no  one  has  actually  analyzed  the  propagation 

9  C 

of  this  instability  in  a  bending,  converging  electron  beam:  all  effort  so  far 

has  been  devoted  to  straight  beams. 

The  theory  developed  so  far  on  the  diocotron  Instability  is  unable  to 

account  for  the  large  noise  content  in  the  long  Kino  gun.  In  addition  to  shot 

noise  near  the  cathode,  a  further  mechanism  of  noise  enhancement  apparently  exists 

in  the  gun.  Some  type  of  feedback  mechanism  is  thought  to  exist,  whereby  diocotron 

amplification  of  the  shot  noise  modulates  the  potential  minimum  near  the  cathode, 

thereby  increasing  the  noise  input  to  the  system.  This  feedback  theory  is  confirmed 

32 

by  the  observation  that  shorter  cathode  lengths  yield  lower  noise  levels  . 


Further  evidence  le  teen  In  a  reduction  In  noise  when  the  cathode  temperature 


is  reduced.  This  "teaperature- United'’  operation  eliminates  the  potential 


win-inniM  in  front  of  the  cathode  surface,  thereby  making  the  feedback 
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mechanism  less  effective  *  .  Similarly,  it  has  been  found  that  placing 


a  grid  between  the  potential  minimum  and  the  bulk  of  the  space  charge  reduces 

S3  34 

the  noise  considerably  *  .  The  grid  prevents  the  electric  field  of  the  noise 


modulation  in  the  space  charge  from  reaching  the  potential  minimum 


The  above  considerations  point  out  a  need  for  an  increased  understanding 
of  the  propagation  of  fluctuations  in  an  electron  beam  with  varying  dc  parameters 
Analysis  of  electron  trajectories  and  the  self-consistent  potentials  would  appear 


to  be  only  part  of  the  solution.  Plasma  instabilities  of  this  nature  must 


necessarily  Involve  a  consideration  of  the  electromagnetic  field  equations.  Given 
an  estimate  of  the  zero-order  quantities,  then,  first-order  perturbations  can  be 


derived  which  will  be  closer  to  actual  fluctuations.  Electromagnetic  solutions,  both 
linear  and  non-linear,  may  be  able  to  encompass  the  computer-analyzed  trajectory 
solutions  and  at  the  same  time  point  the  way  toward  increased  understanding  of 
other  types  of  electron  beams,  such  as  the  hollow  beams  used  in  gyrotrons.  These 
considerations  could  lead  in  turn  to  Improved  design  parameters  not  obtainable  by 
empirical  techniques. 


The  present  study  represents  a  numerical  extrapolation  of  the  work 


of  Knauer  (1966)  and  Buneman,  Levy  and  Linson  (1966)  on  straight  electron  beams 


drifting  in  crossed  electric  and  magnetic  fields.  It  is  anticipated  that  an 


increased  understanding  of  these  simple  cases  will  allow  a  similar  analysis  of 
beams  with  continually  varying  dc  parameters.  We  therefore  begin  with  a  derivation 


of  the  perturbation  equations  for  the  general  case,  then  Introduce  several 


simplifications  for  the  present  analysis 
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4.2  General  Equation* 

The  wave  motions  Co  b«  studied  here  occur  in  electron  besas  driven 
by  crossed  electric  end  magnetic  fields.  In  general,  the  beam  will  be 
considered  "cold",  with  no  thermal  agitation  to  disturb  the  uniform  drift. 


The  treatment  trill  be  non- relativistic,  and  a  quasistatic  approach  will  be 

used.  This  is  permissible  since  one  finds  that  the  contribution  from  the  self- 

v2 

magnetic  field  is  down  by  a  factor  — 2  below  the  contribution  from  the  applied 
magnetic  or  space-charge  fields.  We  therefore  have  the  condition 


7  x  E 


together  with  Poisson's  equation 


It  u 


an  equation  of  continuity 


♦  *  - 


and  an  equation  of  motion 

+  (V  •  7)V  -  -n(E  +  V  x  B)  (4) 

where  n  •  ^  .  In  addition,  we  note  the  definition  of  current  density, 

1  ■  pV  and  we  introduce  the  cyclotron  frequency  o>  ■  nB • 

s  »  c 

The  coordinate  system  to  be  used  is  shown  in  Figure  4.1.  The  general 
case  has  a  beam  esierglng  from  a  given  area  on  a  cathode  into  crossed  magnetic 


and  electric  fields.  The  magnetic  field  is  into  the  paper,  along  the  t  axis,  while 
the  electric  field  is  directed  downward,  along  the-x  axis.  The  beam  then  rises 


WSSBA 


II 

0 

0 

Figure  4.1 

Geometry  of  curving  electron  beam,  shoving  coordinate  systems 
and  orientation  of  electric  and  magnetic  fields.  Note  that 
as  the  beam  bends  it  turn  to  the  +y  direction,  with  the  +x 
surface  on  the  upper  side,  the  ~x  surface  on  the  lover. 

I 

D 

0 

0 
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toward  tha  anode,  converging,  accelerating  and  banding  to  tha  right,  to  the 
+y  direction.  Hanca  baas  velocity,  orientation  and  density  are  all  changing 
In  this  "acceleration"  region.  Tha  baas  than  drifts  uniformly  parallel  to  tha 


y  axis,  which  is  the  case  to  be  analysed  In  this  paper.  In  order  to  sec  up 
the  general  equations,  however,  curvilinear  coordinates  p  along  the  baas  and q 
across  the  bass  are  defined,  such  that  in  the  drift  region  q  becomes  coincident 
with  x,  and  p  with  y. 


The  variables  E,  V,  j  and  p  are  split  into  steady-state  and  time- varying 


parts: 


E  »  E  +  E.  exp[i(fa)t  -  kV  )] 
-  ~o  -i  po 


,  etc. 


Here,  V  is  the  steady-state  velocity  of  the  electron  beam.  The  subscripts  1 
po 

on  perturbed  variables  will  be  understood  in  the  following  equations.  Also,  the 
equations  will  be  written  in  Cartesian  style  although  curvilinear  coordinates 
would  demand  scale  factors  which  are  not  unity.  Such  scale  factors  are  easily 
handled  numerically  but  greatly  Increase  the  analytical  complexity.  The  equations 


will  therefore  be  valid  for  sections  of  the  beam  which  are  small  compared  to 
the  radius  of  curvature  of  the  bean. 

The  following  set  of  five  equations  then  follows  for  the  perturbed 


quantities : 


(“  - k V ' 1  t?) 


V  -  IV 

P  po 


/  dV  \ 

(t?  -  -.) 


+  irtt  (6) 

P 


*•  -  Vi  *  ‘VV  ‘  *-.%  +  lnE,  <7> 

(■  -  kv  •  ‘t?)"  ■  ”p.l5  •  w.%  +  iw(poV  +  i£<boV 


+  ikl  -  0 
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-  ik£  ■ 

P  e„ 


AC  chit  point  It  la  possible  to  solve  the  equations  for  any  one  quantity, 
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such  as  B  or  E  .  or  an  electrostatic  potential  Both  Gould  and  Knauer 
P  9 

31 

solved  for  E  .  However,  Buneman  et  al  found  that  solving  for  V  was  aost 
P  9 

fruitful,  since  various  analytical  difficulties  were  thereby  avoided.  For 
exaaple,  this  makes  Che  equations  iasMdlately  compatible  with  the  boundary 


conditions,  which  are  formulated  most  easily  In  terms  of  V^. 


4 .3  Formulation  for  Straight  Beams 

In  ordar  to  roduca  the  five  equations  to  a  single  one  in  V^,  for  a  straight 
beam,  it  will  be  implicitly  assumed  that  the  perturbation  amplitudes  do  not  grow 
along  the  beam.  Convection  terms,  that  is,  terms  involving  p  derivatives,  will 
therefore  be  dropped,  although  this  does  not  rule  out  the  possibility  of  growth 
rates  appearing  as  complex  k.  Density  0Q  will  be  assumed  constant  and  also  the 
velocity  shear  will  be  taken  constant: 


3V 

T?  ' 


(ID 


Finally,  because  we  are  analyzing  a  beam  in  the  drift  region  of  tube,  subscripts 
p  will  be  changed  to  y,  q  to  x. 

The  equations  then  become: 


(u>  -  kVyo)Vy 

(u>  -  kVyo)Vx 

(u>  -  k^Q)p 
3E_ 


-  ivx(wg  -  cs.)  +  inEy 

-  iu3cvy  + 

•  +  ipoir 


dx 


0 

£ 


Following  Buneman  et  al,  we  substitute  Ey,  and  p  from  the  first  three 
equations  into  the  last  two,  obtaining 


[i(w  -  kVyQ)Vy  +  (u#  -  wc)VxJ 


ik[i(w  -  kV  )V  +  w  V  ] 
1  yo  x  c  yJ 


(12) 

(13) 

(14) 

(15) 

(16) 


(17) 


(1(U)  -  kVyo)Vx  +  a)cVyj'  -  -  kvy0>vy  +  <«,  - 


Is  Che  square  of  Che  plasma  frequency  and  primes  denote 


differendadon  with  respect  Co  x.  The  derivatives  are  taken  inside  the  brackets 


is  eliminated  by  multiplying  the  first  equation  by  w  ,  the  second,  by 


l(u  -  kV  ) ,  and  subtracting.  This  yields 


where 


Substituting  this  expression  for  V  back  into  the  first  equation,  one  finds 


This  equation  may  be  rendered  in  the  form 


where 


Boundary  Condition* 


The  surface  of  Che  beam  is  undulating  because  of  perturbing  surface  waves. 
However,  its  effect  can  be  well  approximated  by  treating  undulations  as  variations 


in  surface  charge  density  of  a  rigid-sided  beam.  This  is  known  as  the  "mean 


surface  approximation".  Conservation  of  charge  flow  in  the  boundary  layer  is 


expressed  by  the  equation 


where  the  upper  sign  refers  to  the  +x  side  of  the  beam,  the  lower  sign,  to  the 


x  side;  o  is  the  surface  charge.  Taking  a  to  vary  as  exp[i(o>t  -  ky)  ]  one  finds 


Across  such  a  surface.  Internal  and  external  fields  match  according  to  the 


relations 


where  subscripts  1  and  2  denote  fields  inside  and  outside  the  beam,  respectively 


Anticipating  that  the  electrical  potential  will  decay  with  distance  from  the 


ik($  e  +$  e  )exp[i(u»t  -  ky)] 


♦  e  ) exp [1 (ut  -  ky)] 


Therefore  the  following  combinations  are  appropriate  to  the  +x  and  -x  sides 


In  order  for  the  solution  (30b)  not  to  exist  on  the  +x  side  away  from  the 


Similarly,  for  the  solution  (30a)  not  to  exist  on  the  -x  side  away  from  the 


Inserting  these  relations  into  the  second  boundary  condition,  equation  (26b) 


and  substituting  for  a 


From  the  first  boundary  condition  (26a),  E 


components  of  the  internal  electric  field  at  the  beam  edge  are  related  by 


Substituting  first  for  E  and  E  from  the  equations  of  motion, (12) and (13) ,  then  for 

x  y 

V  from  (19),  the  boundary  condition  becomes 


where 


4.5  Solution  Scheme 


An  eigenvalue  technique  was  used  to  solve  equation  (22)  together  with  the 


boundary  condition  (34).  This  was  done  in  two  ways,  one  to  find  complex  u  given 


real  k,  the  other  to  find  complex  k  given  real  a).  The  eigenvalue  formulation  for 

ft2 

(i)  versus  real  k  is  produced  by  multiplying  equation  (22)  by  —  2* 


where 


2kV  2(o>  -  kV  ) 
_ s  _ yo 


2kV 

_ s 


Here,  V  has  been  replaced  by  V.  Central  difference  terms  are  substituted  for  the 


derivatives  to  give  an  equation  in  discrete  V ,  j  -  1,...,  N 


(J  ‘  T'’h)Vi  +  f®T  +  *,V)  vj  +  I'")  V  ■ 


Subscript  1  represents  the  -x  edge  of  the  beam,  subscript  N,  the  +x  edge  of  the  beam. 


The  boundary  conditions  (34)  at  these  edges  can  be  formed  by  Taylor  expansion  about 


the  end  points  (Abramowltz  and  Segun  ,  formula  25.3.4  with  p  ■  1  )  to  produce  the 


formulas 


N-2  equations  are  possible  for  N  points  across  the  beam,  with  spacing  h  between  points. 


By  adding  X2V^  and  X  2V^  to  each  side  of  the  respective  boundary  conditions, 
sufficient  equations  are  available  to  form  an  N  x  N  matrix  M,  such  that 


an  initial  value  may  be  chosen  for  X 


value  found  by  the  formulation 


where  eventually  the  difference  between  X.,  and  X.  is  expected  to  be  less  than 


some  small  quantity  e. 

In  practice,  however,  convergence  can  be  slow  by  this  iterative  technique 
It  was  found  more  effective  to  define  a  complex  difference  function 


where  w  is  a  trial  value  of  used  to  produce  X.  ,  ,  which  in  turn  produces 
i-1  •L“A 

X  by  equation  (44),  hence  by  (39).  Zeroes  of  the  function  d^  are  located 
by  various  standard  techniuqes,  the  most  efficient  being  to  bracket  the  aero  by 
searching  over  w  ^  space,  then  to  locate  it  rapidly  by  quadratic  interpolation 
The  eigenvalue  formulation  for  k  vs.  real  oj  is  produced  directly  from 


equation  (22) 
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where  r"  •  r  end 


a”*  ■ 


V  r 

?* 


2(o)  -  kv  )2  1  2k 20)  (o)  -o)  ) 


s _ c  s 


X2  -  k2 


The  difference  equation  is 


(1  -  I  r"h)Vj-i  +  (-2  +  8"2h2)Vj  +  (1  +  I  r"h)V  +1  -  h2X2V 


and  the  boundary  conditions  are  as  before. 


In  the  numerical  work  presented  here,  eleven  equally-spaced  points 
were  used  across  the  beam,  with  nine  on  the  interior  and  two  on  the  edges. 

In  all  cases  a  close  match  was  found  between  numerically  derived  curves  and 
parts  of  previous  analytical  curves.  No  significant  improvement  was  found  when 
the  number  of  points  across  the  beam  was  increased.  This  is  because  the 
fundamental  mode  and  the  first  harmonic,  adequately  described  by  11  points, 
were  the  principal  sources  of  information  for  the  various  curves. 


? 


Beam  with  No  Velocity  Shear 


We  begin  our  analysis  by  examining  a  uniform  electron  beam  of  finite 
width  2a  drifting  in  crossed  electric  and  magnetic  fields.  Following  Knauer 


we  shall  treat  the  case  of  zero  velocity  shear  across  the  beam.  We  begin  with 


equation  (36),  setting  u>  *  0  in  the  terms  r'  and  s'  to  obtain  r'  *  0,  s 


and  hence 


Several  solutions  follow  from  this  equation.  Setting  ft2  *  0,  we  have  the 


dispersion  relation 


According  to  Knauer,  this  relation  comprises  various  oscillations  and  waves  confined 


to  the  interior  of  the  beam.  These  modes  will  not  be  treated  here.  Other  modes 


which  increase  towards  the  surface  inside  the  beam  are  given  by 


Because  equation  (52) 


Solutions  to  this  equation  are  of  the  form  V 


is  a  free  space  field  equation,  the  charge  density  within  the  beam  remains 


undisturbed.  The  resulting  waves  are  then  surface  waves,  and  the  solution  V 


is  appropriate  to  the  +x  edge,  while  the  solution  V  *  eje  is  appropriate  to  the 
-x  edge.  Inserting  these  solutions  into  the  respective  boundary  conditions,  we 


This  simplifies  Co 


w  2  •  *  2(o)  -  kV  )  [w  *  (w  -  kV  )  ] 

p  yo  c  yo 


(54) 


Setting  u  ■  w  -  kV^Q,  solutions  are  given  by 


u,  •  +  -r  w  *  T  +  2w  2 

l  2  c  2  c  p 


(55) 


“2  2  /“c!  +2“/ 


U>„ 


The  corresponding  curves  of  —  versus  ka  are  shown  in  Figure  4. 2  where  — “  0.5, 

"c  ‘“c 

and  the  normalized  velocity  V  /w^a  «  0.5.  Figure  4.3  shows  the  same  curves 

on  log- log  scale  for  comparison  with  later  diagrams.  In  practical  terms,  the 

curves  may  correspond  to  a  magnetic  field  of  114  gauss,  o>c  ■  2  GHz  cyclotron  frequency, 

a  beam  width  of  0.1  am  and  a  beam  velocity  V^q  -  10s  m/sec. 

The  four  curves  were  confirmed  using  the  same  computer  program  used  later 
to  develop  more  complicated  dispersion  curves.  The  curves  of  Figures  4.2  and  4.3 
are  symmetrically  placed  about  the  line  u>  •  kV^Q  (called  the  "synchronous  axis") . 

The  modes  above  the  synchronous  axis  are  propagating  faster  than  the  beam 
velocity,  while  those  below  are  propagating  slower.  Relative  to  the  moving  beam, 
the  curves  closer  to  the  synchronous  axis  may  be  considered  "slow"  modes,  while 
those  farther  away  are  "fast"  modes.  Knauer  describes  the  particle  motion  in 
the  slow  mode  as  circular,  with  electrical  and  magnetic  forces  opposing,  while 
the  fast  mode  consists  of  circular  particle  motion  in  which  electrical  and 
magnetic  forces  are  aligned.  When  a  velocity  slip  is  present.  Interaction  between 
the  slow  waves  of  the  two  surfaces  of  the  beam  gives  rise  to  the  diocotron,  or 
slipping  stream  instability. 

Also  marked  on  the  diagram  are  the  eigenvectors  which  correspond  to  the 
various  curves.  As  described  in  the  solution  to  equation  (52),  surface  waves  for 


Figure  4.2 

Dispersion  curves  for  vsrlous  nodes  of  oscillation  in  a  drifting 
electron  bean  with  no  velocity  shear.  w£/u>£  “  °*5*  Velocity 
distributions  for  the  nodes  are  sketched  near  curves. 


Figure  4. 3 


Dispersion  curves  as  In  Figure  4. 2 on  log-log  scale.  Those  curves 
above  Che  line  u>  •  ltV„0  represent  waves  travelling  faster  then 
the  bean  velocity,  those  be:  im  the  line,  waves  slower  than 
beaa  velocity. 


ka  «  1  arc  functions  which  decrease  exponentially  froa  the  surface  into  the 
been.  For  higher  ka,  various  types  of  Internal  perturbation  are  evident. 


4.  7  Beam  with  Velocity  Shear.  No  Drift 


For  long  wavelength  instabilities,  lea  «  1,  the  flow  nay  be  considered 


incompressible.  Qualitatively,  this  is  related  to  the  tendency  of  particles 


under  such  perturbations  to  behave  as  they  do  in  equilibrium,  and  drift  at 


the-E/B  velocity 


field  perturbations  are  then  divergenceless  as  well 


as  curl-free,  and  so  are  the  velocity  perturbations  (being  the  E  pattern  turned 
through  90°  The  two  boundaries,  however,  undulate  under  the  perturbing  field, 
and  represent  surface-charge  waves  that  communicate  with  each  other  across  the 
apparently  "empty"  space  between  them.  For  a  convenient  description  of  this 
interaction,  one  may  choose  a  frame  of  reference  in  which  the  centre  of  the  beam  is  at 
rest.  Since  this  frame  is  symmetrical  with  respect  to  the  two  interacting  surface 


waves,  the  resultant  growing  wave  mode  can  be  expected  to  be  stationary.  Its 
variation  in  time  and  space  is  then  expressible  as  exp  (lut  -  iky)  where  u)  is  pure 


For  the  case  u  ■  0,  no  imaginary  solution  for  w  is  expected  since 


equations  (55)  do  not  permit  such  a  solution.  Buneman  et  al  has  treated  the 
case  u  f  0  for  the  special  circumstance  u  u>  •  w  2  .  This  relation  follows 


from  the  dc  electrostatic  fields  set  up  by  the  presence  of  the  drifting  electrons 


Such  fields  are  directed  Inward  toward  the  beam,  in  the  x-dlrectlon,  according  to 


where  p  <  0  is  the  charge  density  of  the  electron  beam.  Under  the  influence 


of  this  secondary  B'  field  and  the  ambient  magnetic  field,  the  electrons  drift 


with  a  velocity  -E'/B.  The  velocity  shear  is  then  given  by 


B unman’ s  results  ars  shorn  In  Figure  4.4.  The  velocity  shear  Is 
(D  *  0.5u  in  this  case  (i.s.  6  *  0.5).  The  pure  imaginary  curve  is 
apparent  for  low  ka.  As  ka  Increases  fro®  these  values,  the  two  interacting 
surface  waves  become  further  removed  from  each  other  and  growth  rate  diminishes. 
This  is  a  consequence  of  the  incompressible  nature  of  the  beam  interior:  since 
electric  fields  are  derivable  from  a  potential  that  solves  Laplace's  equation, 
and  since  variations  are  sinusoidal  in  the  y  direction,  the  time  and  space 
variation  must  be  of  the  form 

<(»  *v>  e-**  exp[l(b>t  -  ky)  I  (58) 

This  exponential  behaviour  in  the  x  direction  was  alluded  to  earlier.  Lord 

Rayleigh's  derivation  of  the  ka  «  1  situation  shows  that  above  ka  -  0.64  the 
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instability  becomes  stable  and  that  a  real  solution  exists  (eg.  Hasegawa  ). 

u 

Buneaan's  long  wavelength  derivation  for  general  8  *  ~  <  1  shows  a  transition 

c 

from  instability  to  stability  at  ka  given  by 

ka  -  0.64  +  0.326  +  0.02B2  (59) 

or  ka  ■  0.805  in  our  case  (8  ■  0.5).  The  maximum  growth  rate  for  ka  «  1  is 

given  by 

Im—  -  0.206  +  0.16  -  0.0462  (60) 

wc 

For  intermediate  wavelengths  no  unstable  waves  are  present.  For  short 
wavelengths,  above  ka  given  by 

ka  -  0.5B"X  +  0.25  +  0.1256  -  0.11B2  (61) 

or  ka  -  1.285  in  our  case,  growing  waves  are  again  present.  The  instability 
is  due  to  resonance  at  the  gyrofrequency  between  a  surface  wave  and  a  compressive 


Figure  4. 4 


Dispersion  curves  derived  analytically  by  Buneman,  Levy  and 
Linson^1  for  a  stationary  beam  with  velocity  shear  m9  •  0.5  wc. 
As  In  Figure  4.2|WJ/U  |  •  0.5.  The  real  u>  curve  is  shown  as  a 
solid  line,  the  pure  imaginary  u>  curves,  as  dashed  lines.  Note, 
the  rightmost  imaginary  id  is  a  factor  of  10  less  than  shown  here 


wave  within  the  beam.  At  these  wavelengths,  the  previous  assumption  of  an 

Incompressible  beam  is  no  longer  valid.  The  short  wavelength  instability  occurs 

for  beams  wide  relative  to  the  wavelength  (specifically,  for  which 

k(2a)co  >  to  Is  true) .  The  surface  wave  propagates  without  strong  coupling  to  any 
s  c 

other  wave;  it  sees  small  reflections  from  the  gyro-resonant  layer  through 
a  substantial  width  of  almost  transparent  beam.  The  dispersion  relation  for  this 
mode  is  given  by 

~  -  [kaB  -  0.56  -  0.2562]  +  i(ff)  e"2/B  (62) 

c 

The  real  and  imaginary  parts  are  shown  in  Figure  4.4. 


Numerical  Analysis  of  Beam  with  Velocity  Shear,  No  Drift 


Dispersion  curves  were  developed  numerically,  as  described  in  section 


The  curves  are  shown  in 


Figure  4.5.  They  are  similar  in  several  ways  to  Buneman's  curves.  Figure  4. 4, with 
the  one  striking  difference  that  the  real  curve  of  Figure  4.4  has  been  crossed  by 


a  second  real  curve  running  approximately  linearly  from  u  ■  u>  on  the  u>  axis  to 


in  the  division  of  the  two  curves  into  an  upper  branch  and  a  lower  branch.  Figure 
4.6  shows  the  disperslou  curves  in  a  linear  plot. 


imaginary  w  curves  below  ka  ■  0.77 


are  Identical,  as  are  the  corresponding  real  dispersion  curves  away  from  the 


cross-over  point.  To  the  right  of  the  lower  branch,  a  complex  curve  occurs  in 
Figure 4.  5  which  has  a  corresponding  complex  curve  in  Figure  4.4.  The  continuation 
of  this  latter  curve  is  not  certain,  as  the  points  were  difficult  to  locate  beyond 
ka  -  1.8.  The  imaginary  curve  is  shown,  with  the  real  parts  being  less  by  a 
factor  0.01.  Three  points  on  this  curve  are  given  in  the  following  table: 


(0.000115,  0.05670) 


(0.00112,  0.0900) 


(0.00231,  0.0941) 


In  Buneman's  analysis  an  asymptotic  value  of  0.005  was  given  for  the  imaginary 


at  8  •  0.5,  for  large  ka.  The  present  analysis  suggests  a  larger 


growth  rate  for  gyro-resonance  conditions  in  straight  electron  beams.  This  may 


be  related  to  the  high  gain  observed  in  gyrotrons.  In  addition,  the  real  parts 


are  less  than  in  Buneman's  formula,  equation  (62).  That  equation  implies,  in 


»H» 
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Figure  4.5 


Dispersion  curves  derived  numerically,  using  the  same  parameters 
as  in  Figure  4. 4(statlonary  beam  with  velocity  shear  <o8  -  0.5  <oc) . 
Real  to  curves,  corresponding  to  stable  waves,  as  shown  as  solid 
lines,  while  Imaginary  curves,  corresponding  to  growing  modes 

stationary  in  the  rest  frame  of  the  beam,  are  shown  as  dashed  lines 


Figure  4 . 6 


Dispersion  curves  for  a  saeionary  beam  wlch  velocity 
shear  w  -  8u  .  Curves  are  shown  for  6  ■  0.3.0A  and 
0.5.  Tne  dashed  lines  represent  imaginary  at,  to 
which  the  corresponding  real  part  is  a  fraction  0.01 
or  less.  Note  the  tendency  to  higher  ka  as  6 
decreases. 


-  37 


face,  chac  the  real  part  Is  none  other  than  the  real  branch  of  the  dispersion 
curve.  Figure  4.4.  No  such  complex  solution  was  found  in  our  analysis.  Another 
difference,  a  minor  one,  is  that  the  point  of  commencement  of  our  complex  curve 
is  at  ka  ■  1.45,  rather  than  the  value  ka  •  1.285  given  by  Buneman  et  al 
(equation  (61)). 

The  physical  meaning  of  the  real  dispersion  curves  in  Figures  4.5  and  4.6 

<■  Srf. 

can  be  found  from  an  examination  of  the  associated  eigenvectors.  These  are 
shown  in  Figures  4.7  and  4.8  for  various  points  along  the  curves.  Small  sketches 
of  these  are  also  shown  in  Figure  4.5.  The  upper  branch  is  apparently  a 
surface  wave  mode,  while  the  lower  branch  is  a  bulk  wave  mode.  As  ka  increases, 
the  relative  depth  of  penetration  of  the  surface  wave  into  the  beam  interior 
decreases.  Note  that  in  the  reference  frame  in  which  the  centre  of  the  beam  is 
at  rest  all  waves  are  propagating  in  the  +y  direction,  but  the  electrons  on  the 
-x  side  are  moving  in  the  -y  direction  while  those  on  the  -He  side  are  moving 
in  the  +y  direction. 

The  two  intersecting  dispersion  curves  may  be  related  to  the  "slow"  and 
"fast"  waves  of  Figure  4.3,  where  velocity  shear  is  absent.  Figure  4.6  shows  the 
curves  of  Figure  4.5  on  linear  scale,  for  8  ■  0.5,  0.4  and  0.3.  As  8  decreases, 
it  appears  that  the  dispersion  curve  which  forms  the  left  half  of  the  upper 
branch  and  the  right  half  of  the  lower  branch  is  shifting  to  higher  ka.  In  fact, 
it  appears  to  be  rotating  about  some  point  just  to  the  left  of  the  w  axis.  The 
other  dispersion  curve  is  likewise  turning  toward  the  k  axis.  By  an  extrapolation 
of  this  trend  one  could  imagine  a  mapping  into  the  fast  and  slow  curves, 
respectively,  of  Figure  4.3. 


Figure  4.7 

Velocity  distributions  across  the  beam  for  a 
succession  of  points  along  the  upper  branch  of 
the  real  dispersion  curve.  Figure  4. 5.  As  ka 
Increases,  these  show  a  transition  from  a 
surface  wave  on  the  -x  surface  to  one  of  the 
+x  surface 


Figure  A. 8 


Velocity  distributions  across  the  been  for  a  succession  of  points 
along  the  lower  branch  of  the  real  dispersion  curve.  Figure  4.5. 
These  represent  bulk  conpressional  waves  on  the  beam  interior 
(the  velocity  peaks  shown  here  being  the  amplitude  of  waves 
propagating  along  the  beam  with  wavenumber  k,  velocity  u/k) . 


;  u 

f  EJ 
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The  dependence  of  several  curve  parameters  on  6  can  be  ascertained  from 
Figure  4.6.  Consider  the  changeover  point  from  unstable  to  stable  waves, 
given  by  equation  (59).  The  point  has  been  worked  out  for  three  values  of  6. 

The  following  table  gives  a  comparison  between  previous  analytical  values  and 

values  numerically  derived  in  the  present  analysis: 

ka 

8  analytical  numerical 

0.3  0.738  0.98 

0.4  0.771  0.83 

0.5  0.805  0.76 

We  found  that  the  changeover  point  moved  to  lower  ka  with  increasing  8, 
contrary  to  Buneman's  findings. 

The  maximum  growth  rate  for  ka  «  1  is  given  by  equation  (60) . 

The  following  8  dependences  were  found: 

max  Im  ^/u>c 

8  analytical  numerical 

0.3  0.113  0.063 

0.4  0.117  0.084 

0.5  0.120  0.123 

Our  results,  while  being  more  strongly  dependent  on  8,  show  a  trend  to 
increasing  Im  w/Uc  with  Increasing  8. 

The  onset  of  growing  waves  in  the  gyroresonance  regime  Is  given 
analytically  by  equation  (61).  The  following  8  dependences  were  found: 

ka 


8 

analytical 

numerical 

0.3 

1.944 

3.09 

0.4 

1.532 

2.10 

0.5 

1.285 

1.55 
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Again,  tha  bahavlour  with  6  of  the  two  results  is  similar. 

The  slope  of  the  real  curve  is  given  by  the  first  derivative 
of  equation  (62)  with  respect  to  ka.  This  is  simply  6.  The  following 
values  were  found: 

d(w/u>c) 

•loIie  -JEST 

6  analytical  numerical  for  u>/wc  *0.1 


0.3 

0.300 

0.30 

0.4 

0.400 

0.38 

0.5 

0.500 

0.50 

The  agreement  here  is  very  close.  The  only  strong  disagreement,  then, 
between  previous  analytical  and  present  numerical  results  is  the  dependence 
with  6  of  the  changeover  between  unstable  and  stable  waves  in  the  region 


ka  <  1 
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4.9  Beam  with  Drift  and  Velocity  Shear 


The  cases  presented  in  the  last  two  sections  may  be  combined  to 


produce  a  picture  of  the  behaviour  of  a  drifting  beam  with  velocity  shear. 


Recall  that  the  beam  is  "cold",  with  electrons  drifting  parallel  to  the  beam 


axis.  In  section  4.8  the  beam  centre  was  taken  to  be  stationary.  A  pure 


imaginary  w  was  found  for  ka  «  1;  this  means  that  a  growing  wave  mode  is 


stationary  in  the  rest  frame  of  the  beam.  This  follows  from  the  symmetry  of 


the  oppositely  directed  electron  streams  about  the  beam  centre.  In  a  moving 


beam,  therefore,  one  would  expect  cu  »  kV^o  to  a<*ded  t0  th®  rea^  Part  of  this 


and  all  other  dispersion  curves. 


In  section  4.7  the  beam  was  taken  to  be  uniformly  drifting.  This  resulted 


in  dispersion  curves  which  were  symmetrically  placed  about  u>  ■  kV^Q.  By  analogy. 


one  might  expect  similar  behaviour  when  velocity  shear  is  combined  with  beam  drift. 


Figure  4.9  shows  the  real  parts  of  the  dispersion  curves  for  8  ■  0.5,  and 


—^2.  >  o.5,  derived  numerically  as  in  the  previous  figures.  As  expected,  u  ■  kV 

V 


has  been  added  to  the  real  part  of  all  curves;  also,  a  second  set  of  curves 


appears  which  are  vertically  symmetrical  about  the  line  u>  *  kV  .  *n  c^e  re8ion 


ka  <  0.65  the  two  real  curves  coalesce  to  form  a  curve  along  the  <u  •  kVyQ  axis; 


the  imaginary  part  of  this  curve  is  precisely  as  in  Figure  4.5.  (It  is  coincidental 

V 

that  the  w  ■  kV  curve  is  asymptotic  to  the  Im  o>  curve,  since  both  8  and  —^7 
yo  wca 

are  0.5.)  As  for  the  complex  curves  to  the  right  of  the  lower  branch  in  Figure  4.5, 


the  equivalent  curves  are  believed  to  exist  to  the  right  of  the  real  curves  on  the 


b)  ■  kV  uls  in  Figure 4. 9  but  were  not  found  numerically, 
yo 


The  case  discussed  here  is  an  intermediate  one  between  that  of  sections 


4.7  and  4.8.  The  relation  to  the  stationary  beam  with  velocity  shear  has  been 


discussed.  The  relation  with  Figure  4.3  of  section  4.7  is  more  complex.  The 


"slow”  curves  have  collapsed  to  the  u  ■  kV^Q  axis  while  giving  rise  to  a 


Figure  4.9 


Dispersion  curves  for  s  drifting  been  with  velocity  shear,  with 
Wp/u>g  •  0.5,  Vpo/wca  »  0.5  and  8  •  «8/ u>c  -  0.5.  The  curves  are 
derived  from  Figure  4.5  merely  by  adding  u  ■  kV  to  the  real  part 
of  each  curve;  curves  are  symmetrical  about  the yIine  tu  ■  kVy0 
for  a  given  value  of  k.  Imaginary  curves  are  not  shown. 
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0 

u 

growing  wove.  Near  ka  •  0.7  they  leave  the  synchronous  axis  Co  produce  a 
symmetrical  pair  of  dispersion  curves.  The  "fast"  curves  interact  the  slow 
|  one,  then  proceed  past  lea  •  1.5  to  produce  the  gyroresonance  curves.  Their 

real  parts  likely  remain  close  to  the  synchronous  axis. 
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4.10  Application  to  Binding  Beams 

The  dispersion  curves  presented  here  show  the  behaviour  of  wave 

■odes  In  a  straight  drifting  or  stationary  electron  beam,  with  or  without  velocity 

30 

shear.  They  represent  an  extension  on  the  work  of  Knauer  and  Buneman,  Levy 

31  , 

and  Linson  .  In  particular,  the  requirement  w  u  »  w  has  been  eliminated  in 

9  C  p 

order  to  allow  beams  with  arbitrary  velocity  shear  u>  to  be  analyzed.  As  a 

9 

beam  bends  in  crossed  electric  and  magnetic  fields,  it  undergoes  a  range  of 

velocity  shears  which  in  turn  give  rise  to  a  variety  of  wave  modes.  For  example, 

a  beam  emerging  from  a  cathode  In  the  -x  direction,  and  bending  toward  the 

+x  direction  (Figure  4.1)  starts  with  zero  velocity  shear,  goes  through  a  maximum 

u  ,  then  approaches  a  fixed  w  as  it  enters  the  drift  region  of  the  tube. 

8  ® 

The  trends  in  dispersion  curve  shape  which  were  discussed  in  connection 
with  Figure  4.6  can  be  extrapolated  with  the  aid  of  Figure  4.9  to  describe  the 
progression  of  wave  modes  as  an  electron  beam  undergoes  changes.  Figure  .4.10 
shows  the  changes  schematically.  First,  as  the  electrons  accelerate,  the 
synchronous  axis  rises  from  a  horizontal  position  to  an  angle  with  the  k  axis 
corresponding  to  the  beam  velocity,  w/k  *  V^q.  Secondly,  the  velocity  shear 
causes  the  two  intersecting  dispersion  curves  to  go  from  a  shallow  angle  to  the 
synchronous  axis  ,  to  a  steep  angle,  and  then  to  some  constant  angle.  Beam 
acceleration  will  cause  the  value  of  velocity  shear  to  be  greater  in  regions  of 
higher  velocity,  but  as  the  beam  straightens  out  and  enters  the  drift  region  the 
velocity  shear  becomes  fixed.  Finally,  beam  convergence  means  that  the  beam  width 
will  be  smaller,  hence  k  greater  for  a  given  part  of  the  dispersion  curve.  That 
is,  the  reduced  beam  width  causes  decreased  wavelength. 


The  following  picture  then  emerges  as  the  beam  advances  from  emitter  to 
drift  region.  Consider  a  frequency  u  where  two  surface  waves  exist,  as  marked 


li 

u 
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Changes  in  dispersion  curves  as  beam  accelerates,  bends  and 
converges  In  crossed  electric  and  magnetic  fields.  The 
progression  Is  read  from  the  bottom  up.  Circles  Indicate 
the  polnta  on  the  curves  corresponding  to  a  fixed  frequency 
U).  One  point  In  the  upper  diagram  (the  left  one)  corresponds 
to  a  growing  wave.  Diagram  to  the  right  shows  changes  In  the 
wavelength  as  perturbation  travels  up  the  beam. 
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on  Figure  4.10. As  bean  velocity  increases,  the  wavenumbers  k  of  the  two 
surface  waves  approach  a  common  value.  They  then  transform  into  a  pair  of 
bulk  internal  waves,  at  higher  and  lower  k.  As  velocity  continues  to  Increase, 
the  higher  k  wave  begins  to  propagate  slower  than  the  beam  velocity,  while  the 
lower  k  wave  continues  to  propagate  faster  than  beam  velocity.  As  the  beam 
bends  'into  the  drift  region,  acceleration  decreases .  The 
two  bulk  waves  approach  a  common  velocity  at  the  speed  of  the  beam,  then 
become  a  growing  surface  wave  with  a  large  internal  component.  At  higher 
frequencies  bulk  and  surface  waves  also  exist. 

It  is  recognized  that  a  plot  of  complex  ka  versus  real  u>/u>c  would  be 
Invaluable  in  determining  the  behaviour  of  wave  modes  at  various  frequencies. 

Such  a  plot  was  made  by  fixing  u>  and  varying  k.  All  real  curves  were  confirmed 
as  in  the  previous  figures,  but  no  solution  was  found  with  non-zero  imaginary 
ka.  Such  a  solution  would  not  be  expected  for  a  stationary  beam,  but  could  be 
expected  to  exist  in  a  drifting  beam. 

In  extending  the  picture  just  described,  dispersion  curves  with  constant 
velocity  shear  ui#  can  be  a  function  of  x,  possessing  first  and  second  derivatives 
u£(x)  and  u£(x).  Then  the  equation  in  V  ,  (22)  remains  but  equations  (23) 
and  (35)  take  the  form 


r  ■ 


»2  - 


2kto  (u>  -  kV  )  -  0)  w' 

8  yo _ c  s 

n2 
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Since  kinetic  energy  j  mV2  increases  linearly  with  distance  from  the 
cathode,  velocity  shear  is  approximately  quadratic  across  the  beam.  The 
dispersion  curves  have  not  been  worked  out  for  equations  (63)  and  (64) ,  but  are 
believed  to  be  similar  in  most  respects  to  those  already  derived. 

The  dispersion  curves  presented  in  this  paper  have  been  developed  using 
the  electronic  equations  with  several  simplifications.  A  major  simplification 
has  been  the  removal  of  the  convective  terms,  implying  that  the  perturbation 


waves  maintain  a  constant  amplitude.  If  this  restriction  could  be  eliminated, 
a  quantitative  picture  of  the  development  of  perturbations  in  bending,  converging 
beams  would.be  possible.  Nevertheless,  studies  such  as  the  present  one  are 


necessary  to  establish  the  behaviour  of  wave  modes  in  drifting  electron  beams, 
and  to  make  qualitative  extensions  to  electron  beams  with  varying  dc  conditions 
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5.0  PUBLICATIONS  RESULTING  FROM  THE  AFOSR  SUPPORT 

The  first  year’s  effort  was  Issued  as  a  scientific  report,  MPB-112-1, 
entitled  "Velocity  Noise  in  a  Crossed-Fleld  Diode"  and  was  published  in  Int. 
Journal  of  Electronics,  Vol.  47,  No.  1,  pp.  17-40  (1979).  Appendix  I  contains  a 
copy  of  this  article.  The  second  year's  effort  is  contained  in  Section  4  of 
this  final  report,  MPB-112-2. 

6.0  PERSONNEL 

The  principal  investigator  for  this  study  program  is  Dr.  I.P.  Shkarofsky. 
He  was  assisted  by  Dr.  I.P.W.  Sinclair,  an  industrial  post-doctoral  fellow. 

7.0  INTERACTIONS 

Consultations  on  noise  in  crossed-beam  devices  were  held  with  Prof. 

P.A.  Lindsay  on  September  24,  1978,  of  King's  College,  London,  U.K.  Prof. 

Lindsay  is  one  of  the  original  workers  in  this  field  and  has  written  many 
articles  analyzing  electron  motion  with  a  kinetic  approach.  A  useful  interchange 
of  ideas  arose  out  of  this  meeting. 

Recently,  on  Sept.  24,  1979,  Dr.  G.  Dohler  of  Northrop  wrote  us  a  letter 
expressing  interest  and  asking  for  more  information  on  our  scientific  report 
"Velocity  noise  in  crossed-beam  devices". 

Dr.  V.  Granatstein  of  NRL  has  also  discussed  the  above  work  with  us 


and  a  copy  was  sent  to  him. 
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Velocity  noise  in  a  crossed -field  planar  diode 


I  P.  SHKAROFSKY 


MPB  Technologies  Inc.,  P.0.  Box  160,  Ste.  Anne  de  Bellevue, 
Quebec,  Canada.  H9X  3L5. 


Expressions  are  derived  based  on  kinetic  theory  for  the  seven  moments  over  the 
energy  distribution  function  in  a  planar  thermionic  diode  with  crossed  electric 
and  magnetic  fields,  both  for  the  electric-field  dominated  and  magnetic-field 
dominated  situations.  The  expressions  are  given  as  integrals  over  the  distribution 
function  at  the  cathode  surface,  whose  form  is  left  arbitrary.  For  a  maxwellian 
form,  analytic  results,  valid  in  the  region  of  the  potential  minimum,  are  provided 
in  the  electric-field  dominated  limit  when  the  potential  profile  is  assumed  to  be 
parabolic  in  the  vicinity  of  its  minimum.  The  seven  moments  are  the  density 
profile,  the  two  current  flux  components,  the  two  mean  velocity  components  and 
the  two  mean  square  velocity  components.  From  the  latter  four  moments,  the 
velocity  noise  is  determined.  At  the  potential  minimum,  the  velocity  noise  is 
minimized  when  =  2  where  is  the  plasma  frequency  at  the  potential 

minimum  and  <oc  is  the  cyclotron  frequency,  and  operation  under  this  condition  is 
recommended.  It  is  also  shown  how  Lindsay's,  Ho  and  Van  Duzer’s  and  the 
present  analyses  are  related  to  each  other. 


1.  Introduction 

The  technology  of  crossed-field  devices,  if  understood  better,  could  lead  to 
improved  methods  for  decreasing  their  noise  content,  which  is  one  of  the  main 
problems  at  high  beam  currents.  Although  many  sources  of  noise  have  been 
postulated  in  the  past,  Ho  and  Van  Duzer’s  (1968  b)  theory  is  significant  in  showing 
that  the  random  emission  at  the  cathode,  in  conjunction  with  space-charge  effects 
near  the  potential  minimum,  is  a  capable  initial  noise  growth  mechanism.  Since  their 
analysis  tackled  the  problem  of  current  shot  noise,  we  concentrate  our  attention  here 
on  velocity  noise.  Lindsay  (1960, 1964)  delineates  the  two  regimes  of  magnetic-field 
and  space-charge  or  electric-field  dominated  operation  of  a  crossed-field  diode.  Ho 
and  Van  Duzer’s  (1968  a)  analysis  applies  to  an  electric-field  dominated  case.  Here  as 
in  Lindsay  (1970),  we  discuss  its  relationship  to  Lindsay’s  (1964)  previous  work.  A 
controlled  experiment  performed  by  Bacal  e<  al.  (1970)  shows  the  transition  in  the 
cut-off  characteristic  as  the  current  is  increased  to  an  electric-field  dominated 
condition  from  a  magnetic-field  dominated  condition.  When  it  is  magnetic-field 
dominated  at  low  currents,  the  anode  current  displays  and  abrupt  cut-off  as  the 
magnetic  field  increases  past  a  critical  value,  more  or  less  in  agreement  with 
Lindsay’s  (1960)  theory.  At  high  current  levels,  the  cut-off  characteristic  broadens 
considerably.  Ho  and  Van  Duzer’s  (1968a)  analysis  seems  to  be  more  appropriate 
and  the  less  abrupt  current  decrease  is  consistent  with  Bacal  et  al.'n  results. 

We  employ  a  kinetic  approach  with  an  arbitrary  phase-space  spherically 
symmetric  velocity  distribution  function  at  the  cathode.  Although  numerical  results 
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are  calculated  for  a  maxwellian  distribution,  it  is  just  as  easy  to  do  all  the 
calculations  for  other  distributions  using  the  analysis  given  here.  Implicit  in  the 
generalization  of  the  results  to  an  arbitrary  velocity  distribution  function  is  the 
necessity  to  work  in  polar  coordinates,  in  order  to  reduce,  after  integration  over 
angle,  the  phase-space  integrations  in  three  dimensions  to  only  one  over  the 
magnitude  of  velocity.  For  this  reason,  we  do  not  obtain  error  functions  as  found  for 
a  maxwellian  after  integration  over  two  cartesian  coordinates. 

The  magnetic-field  dominated  situation  is  considered  briefly,  followed  by  a  more 
intensive  analysis  on  the  electric-field  dominated  case.  Various  moments  over  the 
distribution  function  are  derived,  such  as  current  flow  components  and  mean  square 
deviations  in  the  velocity  flow.  In  particular,  we  find  that  the  velocity  noise  summed 
over  both  flow  components  at  the  potential  minimum ,  is  least  when  =  2  where 

cDpai  is  the  plasma  frequency  evaluated  at  the  potential  minimum  and  wK  is  the 
cyclotron  frequency. 


y  CATHODE 


Figure  1.  Geometry  of  a  thermionic  diode  with  crossed  electric  and  magnetic  fields. 


2.  General  analysis 

The  geometry  considered  is  shown  in  Fig.  1.  A  static  electric  field,  E0=  —V<p0, 
modified  by  space-charge  effects,  acts  along  the  x-direction,  and  the  external  field,  B0 
is  along  the  z-direction.  As  in  Lindsay  (1960),  we  define  o)e  =  eB0/m  and  the 
normalized  parameters 

2e<*>0  v  <oex  w#  v 


where  v,  is  a  convenient  normalizing  velocity,  e.g.  vf  =  2lcT/m  for  a  maxwellian 
distribution  at  a  cathode  temperature  T.  Also  e,  m,  and  k  are  respectively  magnitude 
of  electron  charge,  mass  and  Boltzmann’s  constant.  The  relations  between  (wx,  wy) 
within  the  diode  and  their  cathode  values  (irx0,  u>y0)  are 

w^  +  wy~  wlo  -t-  wy0  +  t] 


wy=wy0  +  X  and  wx  =  wx0-X(X  +2u>y0)  +  f/ 

This  basic  starting  point  is  used  in  the  following  sections  to  delineate  the  areas  in 
(Wjto,  m>jo)  space  where  electrons  are  allowed  to  be  present  and  over  which  one 
integrates  the  various  moments.  In  contrast  to  Lindsay  (1960,  1964)  we  map  the 
areas  in  the  (tvx0,  wy0)  plane  rather  than  in  the  (wx,  wy)  plane.  This  is  because  the 
velocity  distribution  function  that  we  use  is  allowed  to  be  of  an  arbitrary  form  which 
is  given  at  the  cathode  position.  We  prefer  to  refer  all  velocities  back  to  the  cathode 
position,  rather  than  carrying  factors  such  as  wj  +  tcj  —  tf  in  the  argument  of  the 
velocity  distribution  function.  As  a  check,  all  the  results  in  the  following  sections 
were  reformulated  and  rederived  in  the  {wx,  wy)  variables,  but  this  analysis  is  omitted 
here. 


Velocity  noise  in  a  crossed- field  planar  diode 


We  uae  the  two-  and  three-dimensional  electron  velocity  distributionf  function 
(evaluated  at  the  cathode),  denoted  respectively  by  / 2(wxo<  wlo)  and  /j(wJo,  «$>, 
wj0).  The  basic  assumption  is  made  that  these  functions  are  spherically  symmetric 
satisfying /j(wJ0,  «$>,  w20)  =/j(«’o)  and  wjo)  =/j(w3-0)  where  w20  =  u>3-0  +  w\0 

and  u>2.0  =  w30  +  w20.  We  define  J2  in  terms  of  f3  as  follows: 

/jfW’xO.  U)Jo)  =  Jt_,/J  j"  /j(«)J0,  w2y0.  wl0)Awrt  (3) 

The  Jt_ 1/2  factor  is  inserted  for  convenience  so  that  for  a  maxwellian  distribution,  f2 
apd  /3  have  identical  forms,  namely 

fi(w o)  =h(n o)  =  exp  ( -  wo)  =  (2 nT/nJ'2  c,3)exp  ( -  w£)  (4) 

where 


vx0dvx0 


The  most  general  moment  at  position  (X,Y)  where  the  potential  is  f]  is  given  by 


ufxdwx  J 

\A  ((WyO  +  ^rdWyQ 


uffdw. 


w’.o  ^.o  {[«>*>  “  X(X  +  2Wyo)  +  »|]“ ~  ‘ }  lllf2(w]_0)  (5) 


using  eqns.  (2)  and  (3),  with  no  variations  considered  in  the  2-direction  of  the 
magnetic  field.  The  square  root  can  be  positive  or  negative  depending  respectively  on 
the  sign  of  *.  For  n=0,  one  obtains  a  wx0/wx  factor  in  the  integrand.  This  is  not 
troublesome  and  no  singularities  appear  after  integration,  as  given  later  in  eqn.  (8). 
The  limits  on  wx0  and  differ  for  the  magnetic-field  and  electric-field  dominated 
cases.  The  following  moments  are  derived,  electron  density  in  configuration  space,  n, 
the  x -  and  y-directed  current  flux  densities,  Jx  and  J,,  the  mean)  velocities,  vx  and  vy 
and  the  mean  square  velocities,  v\  and  v2 .  These  are  defined  by 

n  =  cf3^o,o,  J,  ~  v*  *o,i.  ^  v*/C| =  B02/M 

v^Jvf  =  ffi  0,3/ ^  o,i’  Jy  —  v*B  j,o,  vf/v,=B  2  ,o/^i.o,  (6) 

vr  lvi  =  M  J.oA®  1,0 

The  noise  velocity  components  can  be  shown  to  be  related  to  the  deviations  in 
velocity  by 

<*’*>=2«A/[ti3  ~(vx)2]/(eA,Jx) 


t  Our/,(w)  -  a(r,  w)/»3  =/(w|r)n(r)/t>3  and/3(v)  =  n(r,  v)  =/(v|r)n(r)  where/(v|r)  and  n(r,  v) 
are  commonly  found  in  the  literature,  and  are  referred  to  as  phase-space  density  functions. 

)The  means  vX  f  and  v2  are  actually  related  respectively  to  the  second-  and  third-  order 
moments  over/3(W)  or/(w|r).  They  are  the  first  and  second  moments  over /3(w|r),  also  found  in 
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where  A/  is  the  frequency  bandwidth,  A,  is  the  emitting  area,  and 


are  averaged  over  velocity  and  time.  Equation  (7)  is  based  on  Rack’s  (1938)  analysis, 
a  summary  being  given  in  Smullin  and  Haus  (1959),  and  assumes  that  the  mean 
square  fluctuations  in  the  velocity  components  are  caused  by  shot  noise.  The 
parameters  in  eqn.  (7)  are  to  be  evaluated  at  the  potential  minimum  in  the  presence 
of  space-charge. 

The  double  integrals  in  eqn.  (5)  are  transformed  to  a  single  integral  over  the /2(u>2) 
distribution  function  (omitting  the  i0  subscript  on  to2).  For  this  purposed  we  write 
wx0=w  cos  9  and  wyi)  =  wsinO  and  integrate  over  6,  making  sure  to  cover  the  same 
areas.  The  areas  and  limits  are  given  in  the  following  section.  Here  we  show  the 
results  of  the  9  integrations,  consecutively  for  each  of  the  seven  moments  in  eqn.  (6), 
without  inserting  the  limits.  The  A,(Z)  functions  result  after  integration  over  9  or 
Z  =  sin  9,  with  a=0, 1 ...  ,6,  where 


h„(Z(w ))  =  j  (X  +  wZ)m[w2  +  n-(X  +  wZ)2]( "  - 1)12  dZ 

After  inserting  the  limits,  Z  becomes  a  function  of  w.  We  obtain 

=  [wf2{w2)h„(Z(w))dw 


A0(Z)  =  sin  1  m,  w  =  (wZ  +  X)/(w2 +  ti)112 

ht{Z)  =  wZ,  A2(Z)=i(if-Hc2)[8in_1  m  +  m(1  —  u2)i/2] 

A3(Z)  =  (wZ  +  X)[ri  +  w1  —\(wZ  +  X)2]  -  X(w2  +  q -$X2) 

=h'3(Z)-x(w2+n-\x2) 

A4(Z)=  —  [w2  +  t\  —  (wZ  +  X)2]212 

Aj(Z)=i(f)  +  w2)[sin_,M  — m(1  -m2),/2] 

h6(Z)  =  [tc2  +  >f- (wZ  +  X)2]1/2{|[«c2  +  ri-(wZ  +  X)2]-  w 


In  Aj(Z),  we  use  h'3(Z)  to  designate  the  Z-dependent  part  since  the  other  part  cancels 
after  inserting  the  upper  and  lower  limits. 

The  subscripts,  0, 1 , . . . ,  6  refer  respectively  to  the  moments  n,  Jx,  vx,  v2 ,  Jy ,  vy 
and  v2 .  In  the  following  sections,  the  limits  on  w  and  on  6,  the  latter  giving  Z(w),  are 
provided. 


t  For  simplicity  in  notation,  we  use  as  the  dummy  variable  of  integration  u>  instead  of  w10 
to  designate  (te’£0  +  wro)>/1  ftt  the  cathode.  It  should  not  be  confused  with  its  other  use  to 
designate  the  normalized  velocity  within  the  diode. 
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3.  Magnetic-field  dominated  results 

This  case  is  defined  by  the  condition  that  y>y,  throughout  the  diode,  where 


Subscript  ‘a’  refers  to  values  at  the  anode.  Figure  6  in  Lindsay  ( 1960)  shows  shaded 
areas  of  accessibility  of  electrons  in  a  wy  versus  wx  plot  based  on  eqn.  (2).  We  prefer  to 
work  in  the  {wxQ,  w^)  plane  and  the  equivalent  plotf  is  shown  here  in  Fig.  2,  with 

=-*U  +  Xl  +  2Xtwy0  and  w\ot  =  -ti  +  X2  +  2Xu>y0  (10) 


Figure  2.  Electron  velocity  regions  of  accessibility,  shown  shaded,  on  a  versus  plot, 

when  y>y,,  the  magnetic-field  dominated  case. 

These  represent  the  wx0  values  that  allow  electrons  to  reach  the  anode  at  X,  or  point 
X  respectively.  For  the  moments,  Jx  and  vj ,  n  is  odd  and  the  contributions  from 
regions  3.4  and  4A  cancel  those  from  regions  3 B  and  4 B.  Thus  the  limits  on  the 
integrals  in  eqn.  (5)  are  given  by 


For  the  other  moments,  n,  Jy,  vx,vy,  and  vy  ,  the  whole  area  has  to  be  included.  Rather 
than  cancelling,  the  above  indicated  contributions  add,  giving 


t  Region  1  in  our  Fig.  2  corresponds  to  the  shaded  area  below  u>y  =  in  Lindsay  (1960), 
region  2  to  that  between  wy = wyA  and  wy=wyA,  region  3 A  to  that  bounded  by  wy=wyl,  wyA  and 
wyt  on  the  right-hand  side  and  region  3 B  to  the  same  on  the  left-hand  side,  region  4 A  to  that 
bounded  by  the  wx — 0  axis,  wt = u>.2  and  wyA  on  the  right-hand  side  and  region  4 B  to  the  same 
on  the  left-hand  side.  Areas  to  the  left,  which  image  those  on  the  right,  correspond  to  particles 
turning  around  and  arriving  back  at  the  cathode. 


r 

r 

fr. 

r- 

f.  f 

=  1  dw, 0.. 
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As  outlined  above  in  eqn.  (8),  we  change  to  polar  coordinates  and  cover  the  same 
areas  in  integration  over  w  and  0  with  wx0  =  w  cos  0  and  u>y0  =  w  sin  0.  The  following 
Z(w)  functions  occur: 

Z0A  =  ±  1.  Z2>J(w)  =  |- X  ±(u>  +  ri)ll2]/w 


z6(w ) = [  -  xt + {w2 + n.)  ll2Vw 

and  the  following  positive  definite  limits  on  w  are  needed 


The  detailed  algebra  is  omitted.  We  find  that  eqn.  ( 1 1 )  is  equal  to  (}aD  or  G^E  when 
y,  >0  or  y.  <0  respectively,  where 


w  dwf2(w2  )|  ha(Z0)  -  h,(Z , )  1  + 


w  dwf2(w2)[  ha(Z6)  -  h%(Zx )  | 


Also  the  latter  additional  two  terms  in  eqn.  (12)  are  equal  to 

(i)  Q%'f  —  @%,d  "’hen  y, > 0  or  >lt>X2  and  t\>  X2 

(ii)  GaC  —  when  rj,  <  X2  and  t)>  X2. 

(iii)  OmM-OaE  when  >la<X2  and  - X2 <rf<X2,  and 

(iv)  Gal-GaE  when  t/,<  Xjand  t|<  -X2. 

The  additional  O  functions  are  given  by 


w  dw/2(u-2 )[  2A,(Z2 )  -  A,{Z6 )  -  h,(Z , )] 


wdwJ2(w1)[2h,{Z1)-hx(Zl)-h,(Z0 )] 


wdwf2(w2)[ha[Z0)-A,(Z2)] 


w  dwf2{w 1  )lA,(Z,)-At(Zf)] 


wdwf1(w2)[A,(Z1)-A.(Zi)] 


w  du>f2{w2  )|  A,(Z , )  -  fc„(Z« )) 
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u'du'/z{w2)\ht{Z3)  -  h,(Zt )] 


In  terms  of  the  O  function*,  the  moment*  are  determined  by 


"'"'V  •"ojtomii. 
vl  =  (t’,V n/J,)  0yDlt,  ■/,  = 
vy  =(V1  \/ n/J 

where  D/E  or  F/G/H/l  signifies  the  use  of  one  of  these  subscripts  in  the  appropriate 

case,  and  the  A,  functions  with  a=0 . 6  are  given  in  eqn.  (8). 

The  integrals  usually  require  numerical  integration  for  given  fi(w2)  functions 
such  as  a  max  well  ian.  However,  for  Jx  and  the  maxwellian  distribution  in  eqn.  (4), 
the  following  analytical  result  is  obtained,  being  the  same  both  for  Gt  D  and  Gl  E. 


1  .FlOlHIl 


i.FGIHH 


+  exp(tf.)  1 


where  <l>  is  an  error  function.  This  is  identical  to  Lindsay’s  ( 1960)  formula  and  it 
involves  only  anode  parameters.  As  shown  by  Bacal  el  al.  (1970),  Jx  decreases  with 
B0  abruptly  above  a  critical  magnetic  field  given  by  (Xt  —  rjJXm)>2. 

4.  Electric-field  dominated  results 

The  electric-field  dominated  case  is  characterized  by  y,  >  y  (see  eqn.  (9)).  Lindsay 
(1964)  considers  this  case,  but  he  does  not  make  much  headway  since  his  integrals 
cannot  be  formulated  properly  before  the  space-charge  distribution  is  known.  Ho 
and  Van  Duzer  (1968  a),  assume  a  parabolic  spatial  variation  in  potential  and  derive 
an  expression  for  Jx  based  on  certain  approximations.  Below  we  make  use  of  the 
parabolic  potential  variation  in  the  vicinity  of  the  potential  minimum.  With  the 
same  assumption,  Lindsay  (1970)  provides  results  for  n,  Jx  and  Jy.  Here,  we  obtain 
results  for  all  the  seven  moments  based  on  an  arbitrary  velocity  distribution  function 
at  the  cathode.  In  the  following  section,  we  clarify  the  relationships  between  the 
above  papers. 

Consider  eqn.  (2).  Since  u>;jX),  we  require  wy0^(u>),o)m,x=il^/-^  —X 
+  ("’xo)/-^lmin  The  maximum  value  of  wy0  occurs  at  some  position  A|  which 
depends  on  wx0  and  t}{X).  To  find  this,  equate  dv:y0/dX i  =0  to  yield 

(M&)min=  —  -STf  +  JT  ! )  —  »|(  JT  i ) 

and  (23) 

Ko)m«=*»?'(Xi )  dr,(X)/dX 

valid  for  0<X<Xt.  (For  a  magnetic-field  dominated  case,  A,  does  not  exist.) 
Unfortunately  one  cannot  solve  to  eliminate  Aj  between  (tt^o)***  — “Voe  end 
(*®xo)miii  =  M’jtOC  before  >;(A)  is  known.  The  wx0C  curve  is  the  relationship  between  the 
(a’xoLiit  and  (u>y0)m«  velocities  required  for  electrons  to  turn  around  at  point  A,. 


,4 


24 
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Define  wl0A  and  «’20B  an  in  eqn.  (10)  and  also 

“ioc  =  X  t(JT ,  +  2w,0)  -ri(Xi) 

At  X  =  X  i(«'y0),  we  obtain  by  differentiating  wx0C 
f?wx0\  X{  \  1  W,  .  ) 

(sr5  =— 7V-  +  I5 — -r-H2xl+2wy0-^xl)}\  = 

\CU'yO/X  =  X,  Wxo(X  [ )  (.2u'i0  du'y0  J*=x, 


wxo(Xi) 


(24) 


(25) 


in  view  of  eqn.  (23).  Differentiation  of  eqn.  (10)  for  wx0B  and  wx0A  gives  the  same 
result  at  points  Xt  =  A’  and  X i  =  Xt  respectively.  Sint*  the  slopes  are  the  same,  we 
have  at  least  proven  that  the  ww0C  curve  (i.e.  the  («’}.oL.«  versus  (wx0)min  plot)  is 
tangent  to  the  wx0A  and  wx0B  parabolas  for  any  tf{X )  variations.  As  shown  by  Lindsay 
(1970),  this  is  also  true  in  a  (wx,  wy)  plane,  namely  that  the  curve  is  tangent  to 
both  the  wxA  and  u\B  curves  at  their  meeting  (Miints.  Referring  to  Fig.  4  in  Lindsay 
(1904),  this  points  towards  the  situation  in  his  Fig.  4  (a),  with  w ^  touching,  rather 
than  cutting,  the  wyA  and  wyB  curves.  Also  for  a  parabolic  rj(X )  variation,  the 
curve  reduces  to  a  straight  line. 

As  before,  we  prefer  to  work  in  the  (wx0,  wy0)  plane  rather  than  the  (wx,  wy)  plane. 
The  equivalent  of  Fig.  4  (a)  in  Lindsay  (1964)  is  shown  heref  in  Fig.  3.  Point  3  is 
where  the  wx0C  curve,  drawn  here  as  a  straight  line,  touches  the  wx0B  curve.  Points  7 
and  5  on  the  negative  side  are  where  it  touches  the  wx0A  curve.  For  conditions  in  the 
vicinity  of  the  potential  minimum  where  X  =  Xm,  the  important  integration  regions 
over  velocity  occur  around  point  3.  A  straight  line  portion  is  assumed  here  for  the 
M’,o c  curve,  which  as  we  show  below  is  exact  for  a  parabolic  f/(X)  variation. 

Since  from  eqns.  (10)  and  (24),  point  3  satisfies  the  condition  X  =Xt,the  slope  of 
the  straight  line,  using  eqns.  (23)  and  (25),  is  given  by 

dwy0_wx0(X)  U' _ tf _ I1'2 

dwx0  X  L*  X2  J 

The  coordinates  of  point  3  are  given  in  eqn.  (23),  namely  wx0  =  (Xri’  —  r]  —  X2)'12  and 
wyo~W~ X-  Consequently,  the  straight  line  through  point  3  with  the  above  slope 
has  the  equation 


The  plus  sign  applies  to  positive  wx0  and  the  negative  sign  to  negative  u\0  values. 
Define 


tx=-2q(X)/X  +  ri'(X)  and  fi=[Xti’(X)-'t(X)]/X2  (27) 


t  In  Fig.  4  (a)  in  Lindsay  ( 1964),  draw  an  image  curve  of  j/.y-  on  the  right  extending  from 
wy  =  wy6towy  =  wy}.  Region  1  in  our  Fig.  3  corresponds  to  the  shaded  region  below  tvy  =  wyt  in 
Lindsay,  region  2  A  to  that  between  try  =  irf6  and  wf  =  uiyA,  region  2  B  to  that  to  the  right  of  the 
mv  image  curve  between  wy  =  wy6  and  wy  =  wyS,  region  3  to  that  above  wy  =  wyS  and  to  the  right 
of  the  wyA  curve,  region  4 A  to  that  bounded  by  the  wx  =  0  axis,  the  ut^  image  curve  and  wy 
—  wy i  on  the  right-hand  side  and  region  40  to  its  image  on  the  left-hand  side,  region  4 B  to  that 
bounded  by  the  wx = 0  axis,  wy  =  u>yl,  wy  =  wyi  and  the  image  curve  on  the  right-hand  side 
and  region  4 D  to  its  image  on  the  left-hand  side,  region  5 A  to  that  bounded  by  the  wx = 0  axis, 
wy  =  wy)  and  the  curve  wyA  on  the  right-hand  side  and  region  5 B  to  the  same  on  the  left-hand 
side.  Similar  remarks  apply  to  the  (u>„,  wy)  plot,  Fig.  4,  in  Lindsay  (1970). 


:  /•/  ■%  tv  r*- 

.  *  A.  «  v 


1.  P.  Shkarofalcy 


The  exact  expressions  for  (0  and  £,  are  respectively  =0)  and  —  Xt.  Their 
exact  values  are  shown  later  not  to  be  critical,  since  in  our  calculations,  in  fact,  we 
later  assume  Co-*00  &nd  £,-*oo.  In  the  meantime,  we  proceed  with  the  expressions  in 
eqn.  (31). 

With  reference  to  Fiji;.  3,  and  considering  the  moments  Jx  and  r\ ,  the 
contributions  from  regions 4A,4li and  5 A  respectively  cancel  those  from  regions 4Y ’, 
4/>  and  Mi.  The  limits  on  the  integrals  in  eqn.  (5)  are  then  given  by 


WjrOC 


For  the  other  moments,  w,  Jr  t>x,  vy  and  vy ,  the  contributions  add,  giving  the 
additional  contributions,  equal  to 

fwxor  f*  fwxl)A 


J  C  J  wxO B  J  (a  J  wxOB 

In  eqn,  (8),  we  change  from  cartesian  to  |x>lar  coordinates,  covering  the  same 
areas  with  wx0  =  w  cos  0  and  wy0  =  w  sin  0.  The  following  Z(w)  functions  occur  in 
addition  to  those  in  eqn.  (13): 


The  following  positive  definite  limits  on  m*  occur: 

u,o, *  — »)  =  a X-fiX2  when  t;  >  0 
m>o»  =  «2/40 


ii>04  =  ot2/4  —  aA'(/J  —  1 )  +  /?(/!  —  l  )X 2  j 

wl.  =  ac2/4  -ctXjfl-  !)  +  /?(/?  —  I  )X i  J 

Integration  in  polar  coordinates  of  eqns.  (5),  (32)  and  (33)  yields,  after  lengthy 
algebra,  the  following  result  with  the  ht(Z)  functions  in  eqn.  (8)  for  the  sum  of  the 
three  contributions  in  eqn.  (32),  which  we  denote  by  OxA. 
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The  additional  contributions  in  eqn.  (33)  depend  on  whether  X^Xm'  where  Xm 
=  a/2 P  is  given  by  eqn.  (28)  with  ifj.  * 0  at  the  potential  minimum.  We  find  that  these 
terms  are  equal  to  0tB  —  0a  A  if  X  >  Xm  and  GwC—GtA  if  A-  <  Xm,  where 


wdwJ2(w2)[h,{Z2)-k,(ZA)] 


w  dwf2(w2)[h,{Z4)  -  K(Zi)\  +  2 


wdwf2(w2)[ha(Z4)  —  hm(Z6)]  (37) 


w  dwj2(w2  )[h,(Zs)  -  h%(Z , )]  + 


wdwf2(u>2)[h'(Z2)  -  A.(Zj)]  + 


wdwf2(w2)[h,{Zi)  -  hx(Z4)\ 


w  dwf2(w 2  )[h,(Z3 )  -  h,(Zs )]  + 


wdwf2(w1)[h,(Zi)-h,(Zx )] 


In  terms  of  the  G  functions,  the  moments  are  determined  by 

n  =  vfs/nO0BIC,  Jx  =  ”ty/n0i.A'  Vx  =  (v?s/*/Jx)G2,bic 

v2x  =  (v?yj  n/J  x)G  3. A  Jy=vts/ nG4'Bic  Vy  =  (e,V n/Jy)GiBIC 

»  *  ♦ 

v2y  =(v?y/n/J,)G6'B/c 

where  B/C  signifies  the  use  of  subscript  B  if  X  >  A'm  and  G  if  X  <Xm. 


5.  Relationship  of  our  J,  equation  to  those  of  Lindsay,  and  Ho  and  Van 
Duzer 

For  a  maxwellian  distribution,  the  relation  for  Jx,  obtained  from  eqns.  (4),  (8), 
(13),  (34),  (35),  (36)  and  (39),  reduces  to 


+  exp  (>;.)  1 


where  is  the  error  furction.  This  equation  is  independent  of  X.  We  obtain  the  same 
result  from  Lindsay’s  ( 1964)  eqn.  (7.4),  where  n~ 1/2  in  front  of  the  integral  on  the  last 
line  should  read  2n~ 1,2  and  vice  versa  on  the  previous  line,  namely 

+ -r-  f  /i*  ‘exi,(-(«’,o  +  *,)i  +  'r1]<*M’yO 


+  expto.)[l  -4>toi/2)] 


I  ••  I 
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if  we  can  insert  the  parabolic  form  for  rj  from  equation  (28),  X  t  from  eqn.  (29)  and  let 
ri't  =  ri'(X~0).  Equation  (40)  is  also  given  in  Lindsay  (1970)  upon  noting  that  his 
parameter  o  is  (0  —  1  )112 ,  but  one  has  to  omit  the  ^  typographical  error  in  front  of  the 
middle  term  in  his  eqn.  (4.13).  Our  equation  assumes  a  parabolic  potential  but  allows 
arbitrary  velocity  distribution  functions,  whereas  Lindsay’s  formula  either  does  not 
specify  the  potential  variation  (1960)  or  adopts  a  parabolic  variation  (1970)  but 
always  specifies  a  maxwellian  distribution  function.  If  rj(X)  is  not  given,  X^w^)  is 
unknown  and  the  integral  in  eqn.  (41)  cannot  be  manipulated.  Under  both 
assumptions  of  a  parabolic  potential  and  a  maxwellian  distribution,  eqns.  (40)  and 
(41)  agree.  Ho  and  Van  Duzer  (1968  a)  use  both  assumptions  and  in  addition  they 
address  the  limit  when  the  anode  is  faraway  ( X-*co )  and  need  not  be  considered  in 
the  region  of  the  potential  minimum.  In  this  limit  and  provided  is  somewhat 
greater  than  one,  the  two  terms  in  eqn.  (40)  containing  X,  approach  one,  yielding 


This  is  identical  to  Ho  and  Van  Duzer’s  result  who  use  the  notation  fl  =  o)c(/J  — 1)1,2> 
tfm=  —a2/4fi,  and  2toKxm(m/2kT)111  =  <x/p.  Thus  our  result  is  the  same  as  Ho  and  Van 
Duzer’s  in  the  limit  of  and  with  a  maxwellian  distribution.  Finally,  Lindsay’s 

result  reduces  to  Ho  and  Van  Duzer’s  in  the  limits  of  X,-*oo  and  a  parabolic  spatial 
potential. 

Ho  and  V an  Duzer  ( 1 968  a)  proceed  to  apply  a  further  simplification ,  namely  that 
the  potential  minimum  is  deep  and  close  to  the  cathode,  i.e.  xm-»0.  This  essentially 
means  that  ot2/4  >  1  and  /?>  1  but  a2/4/?  can  be  of  order  one,  since  from  eqns.  ( 1 )  and 
(27): 


«- -2i)JXm=*  -4e<t>J(mxmwev0)  and  /?  =  ~nJXi  =  -2e<fcJ(mxlfal)  (43) 


Equation  (42)  reduces  to  the  following,  which  Ho  and  Van  Duzer  show  to  be  a  good 
approximation: 


Generally  we  require  that  fi  ^  1  -5  and  a  >  3  2  for  the  error  functions  to  approach  one . 
This  decrease  of  Jx  with  magnetic  field  for  an  electric-field  dominated  case  is  much 
slower  than  for  the  magnetic  field  dominated  case,  as  is  indeed  observed  by  Bacal  ei 
al.  (1970). 


6.  Simplifications  when  /)>  1*5,  a >3*2  and  X.»  1  for  the  electric-field 
dominated  case  and  calculated  results 

We  simplify  our  relations  in  eqn.  (39)  in  the  limits  of  small  xm  and  large  x„ 
adopted  by  Ho  and  Van  Duzer  (1968a),  equivalent  to  1-5, 0 3  2  and  X,»l.  In 
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eqn.  (35),  we  retain  wob  and  wM  (which  is  equal  to  tcM  at  X  =  JfB),  but  we  let  Wg, 
-♦oo  and  u?o,-*oo  (equivalent  to  £0-»  —  oo  and  £,-*cc).  Thun  0mA  reduces  to  the  first 
term  in  eqn.  (36),  0l  B  to  the  first  two  terms  in  eqn.  (37)  and  G'fcC  to  the  first  three 
terms  in  eqn.  (38). 

These  expressions  can  be  further  drastically  reduced.  (See  the  Appendix.) 


Define  the  function 


dwwf^w1  —a){w2  —  a)*/2 


All  the  moments  in  eqn.  (39)  can  be  written  in  terms  of  pa  after  inserting  the 
equations  in  the  Appendix  and  noting  that  Wgb  =  —  r/m.  The  following  expressions  are 
obtained: 


n(X>Xj=  -n(X<Xm)  +  4nvfPl  (»/,»))  =  2n vf[px(t\,  qj  -  (t\  -  nm) 1/2  Po(<>,  7»)  I 
Jx  =  -  1  )/P]l,2P2(nM,  *lm) 


v,(X>XJ--v,{X<Xm)+ 


These  are  our  final  simplified  relations,  valid  for  an  arbitrary  fb(wg)  at  the  cathode. 
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When  /j  is  a  maxwellian  distribution,  given  by  eqn.  (4),  we  find  that 


n/»r  -  exp  (ij)[  1  ?♦((»»  -  Jl/2)]  for  X^X, 


(exp  (»»  —  »|-) 


[«xp  (n-nm))Pll2l^^{(V-nm)ll2))±(fl-*1m)ll2/Pil2  for 


v^/vf  =  (t\  -rim)/P+i-  (2fi) 
where 

v2  =  2kT/m 


This  set  of  relations  for  a  maxwellian  distribution  can  also  be  derived  directly  upon 
integration  in  the  cartesian  frame  (wx0,  w^).  The  use  of  the  polar  frame  allows  us  to 
derive  the  relations  in  eqn.  (46)  for  an  arbitrary  /3  distribution  function. 

The  following  observations  are  noteworthy,  subject  to  the  assumptions  in  this 
section. 

(1)  The  expressions  for  J„  and  Jy  are  independent  of  X  or  ij(X)  and  in  fact 

U‘/2,  (J2x+J2)ll2=J.exp(r,m)  (48) 

the  latter  being  identical  to  the  expression  in  the  zero  magnetic  field  case  (omitting 
the  t\m  dependence  on  magnetic  held).  The  relation  for  Jy  in  eqn.  (47)  is  also  given  by 
Lindsay  (1970)  in  his  eqn.  5.5.  His  Fig.  9  also  shows  that  J y  is  independent  of  X  over 
a  large  range  of  X  values,  but  because  he  adopts  /?=a  =  2  with  a  not  large  enough  to 
satisfy  our  conditions,  Jy  decreases  near  X  =0.  When  a >  3-2,  Jy  is  independent  of  X 
even  for  X  approaching,  zero,  provided  the  conditions,  given  at  the  head  of  this 
section  are  also  obeyed. 

(2)  A  single  relation  applies  for  ail  X  both  for  v2  and  v2 . 

(3)  The  relation  for  »  is  identical  to  that  in  a  zero  magnetic  field  case.  This 
expression  can  be  thought  of  as  the  ‘ second  approximation'  to  the  density  profile, 
since  inserting  a  parabolic  potential  (d2ij/dX2  =  2fi)  into  Poisson’s  equation,  d2>i/dX2 
=  (2 e/nuoi)d2<t>0/dx2  =  2e2n/mw2e0,  implies  that  n  is  a  constant  in  the  ‘first  approxim¬ 
ation’  given  by 

n=Pma)2eQ/e2  or  (o2/a>2=P  (49) 

where  (Op  =  (ne2/e0m)>/2  is  the  plasma  frequency!-  Using  the  second  approximation 


f  We  note  that  when  n  —  nT  (see  eqn.  (4)),  p  is  equal  to  k~2»(2/X2,  where  k  and  {  are 
defined  in  Lindsay  and  Goodwell  (1965). 
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Figure  4.  Normalized  electron  density  versus  X/Xm  when  <pm=  —  0-155  V  and  T=  llflOK.or 
»|m=*  -1-66.  The  arrows  indicate  ‘effective’  values  of  P  when  for  example  eJ ,  —  1  69 
x  lO'A/m2,  B0  =  660G  and  r.»9x  l<r6m. 


Figure  5.  The  current  density  flaxes  JJ[J,  exp  (n,)]  and  JJ\J,  exp  (»>,)]  and  their  r.m.s.  sum 
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for  n  and  Poisson’s  equation,  one  can  solve  for  X(tj)  and  thereby  effect  an  iteration 
process  (not  done  here).  Probably  one  would  find  a  magnetic  field  dependence  for  n  in 
higher-order  approximations.  Presently,  we  match  /?  to  ai^,  its  value  at  the 
potential  minimum  point,  with  eqn.  (49).  Our  results  are  qualified  to  apply  only  in 
the  region  of  the  potential  minimum. 

Calculated  results  for  the  seven  moments  in  eqn.  (47 )  are  shown  in  Figs.  4-9  up  to 
X/Xm  values  just  beyond  2,  and  for  the  given  values  <j>m  =  —0155  V  and  T  =  1 160  K, 
equivalent  to  t\m=  — 1-55.  If  for  example,  £0  =  650G  and  x»=9  x  10'*  m,  then  Xm 
=0-549,  a  =  —2>tn/Xm  =  5-65  and  0=  —  rjm/-X^  =  5-15.  For  purposes  of  illustration,  we 
allow  in  some  of  the  figures  ft  values  between  1-3  and  200  and  X  values  down  to  zero, 
although  some  regions  may  fall  outside  of  our  range  of  applicability.  Figure  4  shows 
that  the  electron  density  drops  rapidly  away  from  the  cathode  but  decreases  more 
slowly  beyond  Xm.  Lindsay  (1970)  also  finds  a  similar  behaviour. 

Figure  4  also  shows  that  because  n  is  decreasing,  the  ‘effective  /?’  given  by  eqn. 
(49)  is  also  decreasing.  It  seems  that  for  large  X/Xm»l,  p  can  change  from  P>  1  to 
P  <  1  or  from  an  electric-field  to  a  magnetic-field  dominated  situation  within  the 
same  diode.  An  analysis  allowing  for  such  a  change-over  is  non-existent  at  present. 
Since  Jx  and  Jy  are  independent  of  X,  they  are  plotted  versus  /?  in  Fig.  5.  For  the 
magnetic-field  dominated  case,  or  for  the  electric-field  case  at  low  a  and  P,  Jy  varies 
with  X  as  shown  in  Lindsay  (1962),  Lindsay  and  Goodwell  (1965)  and  Lindsay 
(1970).  The  electric-field  dominated  case  with  03-2  and  P^  1-5  seems  to  differ 
greatly.  Figures  7  and  9  show  that  vx  and  vy  are  symmetric  about  the  Xm  point  and 
vary  by  a  much  lesser  degree  than  vx  and  vy .  Figure  8  indicates  that  vy  becomes  very 
large  as  p  increases,  in  contrast  to  vx  in  Fig.  6  which  decreases  with  p. 

In  Figs.  10  and  1 1 ,  we  plot  the  profiles  of  and  (vy)  defined  in  eqn.  (7).  The 
two  behave  differently  for  small  X  and  large  p.  The  former  becomes  less  negative 
whereas  the  latter  becomes  exceedingly  negative.  Negative  values  indicate  that  vx 
>{vl)in,  which  is  possible  since  the  integrand  over  the  (vx—vx)2  average  contains 
the  product  vjit  which  becomes  negative  for  negative  vx  values.  Similar  remarks 
apply  to  the  vy  averages. 

At  X  =  Xm,  curve  (C)  passes  through  zero  in  Fig.  10  and  curve  (E)  in  Fig.  1 1 ,  their 
P  values  being  chosen  for  this  to  occur.  Figure  12  gives  the  sum  of  the  plots  in  Figs.  10 
and  1 1 .  Curves  (C)  and  (E)  coincide  here.  As  P  decreases,  the  curves  move  towards  the 
left  until  p  is  equal  to  2,  and  then  they  move  back  towards  the  right.  The/? =2  curve 
(D)  gives  the  least  negative  value  of  <c2>. 


7.  Velocity  noise  at  the  potential  minimum 

In  the  electric-field  dominated  case,  the  velocity  noise  at  the  potential  minimum 
is  readily  obtainable  from  eqns.  (46)  for  a  general  velocity  distribution  function  and 
(47)  for  a  maxwellian  distribution.  Combining  these  relations  with  eqn.  (7)  yields 
with  arguments  (qm,tlm)  for  the  p  functions  and  with 
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0  0.8  0.8  1.2  1.6  2.0 


x/x„ 

Figure  12.  Normalized  velocity  noise  sum.  =  +  <»*>,  versus  X/Xm  when  r/m  = 

—  1-55.  Curves  correspond  to  0  values  (B)  (I.  ((’)  (D)  2,  (E)  1-480  and  (F)  1-3. 


eA /  [/  0  V'T/20+lW  16/  0  yil 

2eA'itJ,Pi  (\0  —  1  /  \\  0  )pi  6  \0  -  1  /Pi  J 

+  v/fl<3/>-  l)pJPp2  ~  l«^pl/9pilj 


Omitting  the  slight  dependence  of  the  p  factors  on  P,  the  quantity  within  the 
parenthesis  in  |<c2>|  has  a  minimum  with  respect  to  0  at  0  =  2,  and  this  occurs  when 
<«J)=  <t^>.This  is  true  for  any  velocity  distribution  function.  Then 


For  a  maxwellian 


At  best  this  is  4-23  times  as  large  as  (4  -  Jt),  the  value  for  a  zero  magnetic-field  case.  In 
practice,  a  value  of  0=2  is  recommended,  or  0  at  least  within  the  limits  1-6  to  3. 
Inserting  0  =  2  into  eqn.  (40)  gives 

aijm/w?  =  2  (66) 

This  is  the  desired  relation  for  minimum  velocity  noise  at  the  potential  minimum. 
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In  terms  of  the  alternate  relationships,  P  —  - *lJXm 2  and  a=  the 

conditions  /J=2  and  a. ^3-2  imply  that  we  want  X.^O-8  and  \  rim\  =  2Xm2>  1-28. 

Van  Duzer  (1963)  uses  the  magnetic  field  independent  factors  (4  —  7t)  and  2 
respectively  for  <rx2>  and  .  In  contrast,  our  analysis  provides  the  magnetic  field 
dependent  factors  in  eqns.  (50)-(52). 

Figure  13  shows  the  normalized  variations  of  <ex2>,  <t>,2>  and  <e2)  versus  fl  at 
the  potential  minimum.  The  <»x2>  quantity,  normalized  as  shown,  increases  from 
large  negative  values  and  tends  towards  (4  — ji)  as  P  becomes  large.  The  {vy2} 
normalized  quantity  decreases  to  large  negative  values  as  P  becomes  large.  The  sum 
<»2>  is  least  in  magnitude  at  p  =  2,  decreasing  to  large  negative  values  for  larger  and 
smaller  p.  Since  <r2>  is  associated  directly  with  velocity  noise  at  the  potential 
minimum,  the  main  conclusion  of  this  analysis  is  to  operate  the  diode  by  adjusting  B 
and  T  so  that  <Op*/cj2  =  2  is  satisfied. 


Figure  13.  Variations  of  normalized  velocity  noise  components,  (v2),  (v2)  and  their 
sum  (v2)  versus  ft  at  the  potential  minimum. 

8.  Conclusions 

We  analyse  both  the  magnetic-field  and  electric-field  dominated  cases  in  a 
thermionic  diode  with  crossed  electric  and  magnetic  fields.  Assuming  a  parabolic 
spatial  variation  in  potential  about  the  space-charge  minimum,  we  clarify  several 
aspects  and  obtain  simple  but  important  results  for  the  noise  in  the  velocity 
components. 

In  SS  3  and  4,  we  illuminate  Lindsay’s  (1960, 1964)  considerations  on  the  velocity 
regions  over  which  one  has  to  integrate  the  velocity  distribution  functions  to  derive 
the  moments  integrated  over  the  distribution.  We  show  the  following: 

(a)  The  a>xOC  curve  is  tangent  to  both  the  wx0A  and  u>x0B  curves  in  a  (wx0,  w^)  plot 
(see  eqns.  (10)  and  (24))  under  electric-field  dominated  conditions. 
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(b)  The  regions  and  limits  of  integration  are  mapped  and  theirequivalences  in  the 
(10,,  mi,)  and  («0xO>  i0,o )  planes  are  identified,  both  in  the  electric-field  and  magnetic- 
field  dominated  situations. 

For  the  electric-field  dominated  case,  Lindsay  (1964)  could  not  proceed  before  the 
potential  distribution  was  known.  Upon  assuming  a  parabolic  variation  about  the 
potential  minimum,  we  illustrate  in  §5,  the  following,  also  noted  by  Lindsay  (1970): 

(a)  Our  formula  reduces  to  Lindsay's  (1964),  if  a  parabolic  potential  variation  is 
inserted  into  Lindsay's  equation. 

(b)  The  result  of  Ho  and  Van  Duzer  (1968  a),  also  for  a  parabolic  variation,  agrees 
with  our  result  if  we  allow  the  anode  to  recede  to  infinity,  as  they  do.  Our  formula  is 
more  general  in  that  it  allows  the  velocity  distribution  function  at  the  cathode  to  be 
an  arbitrary  function  of  wXQ  +  wy\  +  wt\. 

General  relations  are  derived  for  seven  moments  in  terms  of  integrals  over  the 
arbitrary  two-dimensional  electron  velocity  distribution  function,  both  for  the 
magnetic-field  and  electric-field  dominated  cases.  For  the  latter  case,  we  effect  in  §6 
major  simplifications  upon  assuming,  as  Ho  and  Van  Duzer  (1968  a)  do,  that  the 
potential  minimum  is  sufficiently  deep  and  close  to  the  cathode  but  very  far  from  the 
anode.  These  conditions  are  satisfied  for  high  current  emission  densities  in  high 
power  tubes.  Then  the  moments  are  expressed  in  terms  of  simple  integrals  over  the 
three-dimensional  distribution  function,  which  are  integrable  for  a  maxwellian  form. 
Profiles  for  the  moments  are  plotted  about  the  potential  minimum. 

The  cut-off  characteristic  of  Jx  versus  magnetic  field  differs  markedly  (see  eqns. 
(22)  and  (44))  between  the  magnetic-field  dominated  case,  where  it  decreases 
abruptly,  and  the  electric -field  dominated  case,  where  it  is  broader  and  indicative  of 
such  a  case  rather  than  being  due  to  a  high  electron  temperature.  These  cut-off 
concepts  agree  with  the  experimental  work  by  Bacal  et  al.  (1970).  We  denote  the 
important  space-charge  parameter  by  P ,  equal  to  —  riJX2  generally,  and  given  by 
for  a  parabolic  potential.  Electric-  and  magnetic-field  dominated  conditions 
are  characterized  respectively  by  P  >  1  and  P  <  1 .  When  our  conditions  are  satisfied , 
namely  P^  1-5,  a  =  —  2t\m/Xm^3-2  and  the  anode  position  a:,-»oo,  than  J x  decreases 
slowly  as  «/a[(^  —  1  )/^))1/2 exp (r/m)  as  the  magnetic  field  (Bccp~112)  is  increased. 

The  new  relation  for  the  J,  current  flux  subject  to  the  same  conditions  on  a,  P  and 
xm  is  derived  in  §6  for  the  electric-field  dominated  case  and  gives  J  y  =  J*P~ 1,2 
exp  (p/m),  independent  of  X.  The  r.m.s.  sum,  (Jx2  +  J2y)lll  =  Jx exp  (r/m)  is  identical  to 
that  in  the  absence  of  a  magnetic  field. 

From  the  relations  for  the  velocity  moments,  vxy  and  vx2,  we  obtain,  using 
Rack’s  (1938)  method,  the  noise  velocity  components  <*£,)  in  terms  of  the  mean 
square  deviations,  vlty—(vxy)2,  and  in  particular  their  values  at  the  potential 
minimum.  The  results  are  illustrated  in  Fig.  13.  In  the  absence  of  a  magnetic  field, 
<tr2>  has  the  multiplicative  factor  (4  —  n).  In  crossed  electric  and  magnetic  fields, 
both  (vx  y)  components  depend  on  magnetic  field  through  the  P  factor.  The  sum 
|<»2>  +  0;>|  has  a  minimum  at  p  =  2  and  its  multiplicative  factor  is  |2N/2(5  — 2n)| 
which  is  4-23  times  worse  than  the  zero-magnetic  field  case.  For  optimum  operation, 
one  should  adjust  P  =  t ai^Jw 2  =  2.  In  all  cases,  the  velocity  noise  is  also  proportional 
to  the  cathode  temperature.  If  the  tube  can  be  operated  with  no  cathode  heater,  this 
broad-band  noise  will  be  greatly  reduced.  This  probably  explains  the  experiments  of 
Bespalov  and  Usychenko  (1976)  and  Bespalov  et  al.  (1976). 
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Appendix 

Simplification  of  expressions  for  case  in  §6. 

We  keep  the  first  term  of  G,  A  in  eqn.  (36).  the  first  two  terms  of  6'„  „  in  (37 )  and  the 
first  three  terms  of  GaC  in  (38).  We  wish  to  transform  these  three  approximated 
relations  to  yield  the  results  in  eqn.  (46). 

Define 

X0  =  X-Xm  and  s2  =  w2~ wlb  =  w2  - a2/4/J  =  w2  +  t]m 
The  following  equality  is  particularly  useful: 

[w2  +  ri-(wZ4  j  +  X)2)ll2=P~'l2s  +  (P-\)ll2\X0\Sign(X0)  if  w2'^w\i 
=  |Ao|(0-l)l/2  +  0'1/2sSign(Ao)  if 

where  Sign  stands  for  the  sign,  so  that  Sign  (A—  Xm)=  ±  1  depending  on  whether 
A<Am.  Each  moment  has  to  be  considered  separately  and  integrated  by  parts. 

It  is  convenient  in  the  intermediate  algebra  to  define  F(w2)  and  FD(x2)  functions, 
given  by  dF(w2)/dw=  —2 wf2(u>2)  and  dFD(w2)/du>=  —2 wF(w2),  to  be  used  in  the 
integration  by  parts.  After  considerable  algebra,  one  obtains  the  following  results  for 
the  approximated  forms  of  eqns.  (36)  to  (38). 


,  =  y/0\X  0|J*« 


dsF(s2  +  w20b)/(s2  +  pX20),  ifXo>0 


G0.c=nF(w20a)-G0'B,  if  Xo<0 


GiM=^y-y2jy*°2f2(* 


2  +  u>ob2) 


°i,B=iP112 


ds  F(s2  +  u>ob)  + 


,,x,dsFD(s  +  wlb) 
o  *2  +  PX20 


+  \X0\(P~2)P  1/2  j  ds82f2(a2  +  wlb),  if  Xo>0 
Gi.c~^FD(w20t)-G1B,  if  A0<0 

°3A = \ {nr)  £ d*  *^2(**'+ + “ 1 )] 

0s.B  =  G*.c  =  2P~i/2  j*  dss2f2(s2  +  wlb) 


f  ®  f 

Oj..  =  (l/2)^'/2W  d8F[82  +  w20b)+  - § 

. Jo  Jo  8 


dsFD(a2  +  wlbf 


+PX20 
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-\X0\(P-2)P'112  j*  dss2f2(s2  +  u>lb),  if  X0>0 


Q$,C - 2 ^ °(W0a)  - ifX<)<0 


06.B = 06.c  =  2p - 1,2 1  ds 8 2f2(s2  +  »2»)| Xl  +  S2(3p  - 1  )/3 P) 

All  the  above  integrals  over f2,  F  and  FD  can  be  transformed  to  integrals  over f3, 
the  standard  three-dimensional  distribution  function  in  eqn.  (3).  The  following 
identities  which  can  be  easily  proven  can  be  substituted  into  the  above  relations  to 
yield  the  results  in  eqn.  (46). 


i; 


dss2/2(s2  +  w^b)  = 


Fd(wL)  = 


#1 

'0D 

du>  w3J3(w2  4-  WQh) 

0 

3>/7t 

8 

j*  dww5f3(w2  +  Wob) 

00 

dww3f3(w2  +  wlb) 

0 

2jn 

• 

j"  dwwf3(w2  +  iCoi>)jj 

1 - 1 

7 

a 

8 

3>A . 

j”  dww*f3(w2 +wla) 

4  Jn 
3 

t 

(™  +  PX2o)312  ov2 
(PX20)1'2  PA° 
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